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The primary purpose of the present series of papers is to 
give more explicit mathematical form to the theory of 
Mulliken and Hund, and in some instances to the alterna- 
tive theories of Pauling and Slater and of Heitler and 
Rumer, for valence in carbon compounds. A critical 
comparison is given of the Slater-Pauling and Hund- 
Mulliken concepts of CH,, which are based respectively on 
localized bonds (“electron pairs’’) with the Heitler-London 
method, and on one-electron wave functions (Mulliken’s 
“orbitals”) for a self-consistent field with tetrahedral 
symmetry. The H-M procedure avoids hybridization of 
the carbon 2s and 2 wave functions, but allows two 
(though never three or more) electrons to accumulate on 
one H atom, as well as up to eight on a C atom. Inadequate 
cognizance is thus taken of the tendency of inter-electronic 
Coulomb forces to keep two electrons apart. The S-P 
procedure avoids this excessive accumulation, but at the 
expense of not letting an individual wave function be of a 


symmetry type (irreducible group representation) ap- 
propriate to a tetrahedral field. In particular, its s-p 
hybridization ‘‘undiagonalizes” the internal energy of the 
C atom. These points are illustrated by explicit exhibition 
of the secular determinant, which is the main new feature. 
Because inter-electronic repulsions make the dynamical 
problem more complicated than a one-electron one, the 
tetrahedral symmetry need only be preserved in the 
properties of the total wave function of the entire system 
rather than that of one electron, but the departures from 
individual tetrahedral symmetry should be less than in the 
Slater-Pauling theory if the Hartree self-consistent field is 
really a good approximation. Thus both the H-M and 
S-P methods, though qualitatively exceedingly illumi- 
nating, have their own characteristic drawbacks from a 
quantitative standpoint unless refined by inclusion of 
higher approximations which ultimately merge the two 
methods but which practically are very difficult to make. 





WO approaches have been opened to the 
subject of directed or orbital valence? in 
polyatomic molecules. One method is that de- 
veloped by Slater? and by Pauling.* Their model 





‘A preliminary account was given at the Washington 
meeting of the American Physical Society, April, 1932, 
(abstract in Phys. Rev. 40, 1037 (1932)). 

* Since we are here confining our attention to theories of 
directed or orbital valence, we do not need to include in the 
present discussion the Heitler-Rumer treatment of carbon 
compounds, in which the directional properties exist only 
In virtue of the repulsions between the atoms attached by 
the central carbon atom. This Heitler-Rumer method must 
not, however, be overlooked and will be compared with the 
other methods in Part III. 

*J. C. Slater, Phys. Rev. 37, 481; 38, 1109 (1931). 

‘L. Pauling, J. Am. Chem. Soc. 53, 1367, 3225 (1931). 


is based on what Hund calls a “‘localization”’ of 
bonds whereby particular wave functions are 
assigned to “‘overlap’”’ particular atoms. To ex- 
plain tetrahedral carbon compounds, the so- 
called ‘‘sp*”” valences and wave functions are 
introduced. Namely, Pauling and Slater show 
that by taking linear combinations of the usual 
2s and 2 wave functions it is possible to build 
four orthogonal wave functions which project in 
space like the four corners of a tetrahedron. The 
Slater-Pauling model lends itself particularly to 
calculations by means of the Heitler-London 
method, in a generalized form given by Slater® 
and by Born.’ 


5M. Born, Zeits. f. Physik 64, 729; 65, 718 (1930). 
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The other approach is that due to Hund® and 
to Mulliken,’ and is a generalization of their 
diatomic configuration theory. In Mulliken’s 
recent model, an electron is not assigned to a 
particular valence but is allowed to move in a 
field of the same symmetry as that of the 
molecule in question; e.g., a field of tetrahedral 
symmetry in CH,. Great simplicity rests in the 
fact that with this model the dynamics are that 
of a ‘“‘one electron system,”’ for instance the two 
center problem in diatomic molecules, the five 
center one in CH,. Electrons are, so to speak, fed 
to the molecule one at a time, and by the 
Pauli principle not more than two of them can be 
digested in the same orbital state. Inter-electronic 
repulsions can, in principle, be included to the 
extent to which they may be represented by a 
Hartree self-consistent field of the given sym- 
metry. If the molecule is symmetrical, several 
solutions of the ‘‘one-electron’’ wave equation 
may have the same energy. The degree of 
degeneracy can be determined from that famous 
paper of Bethe* in which the different repre- 
sentations of the rotation group are resolved into 
their irreducible portions in various crystal- 
lographic groups. In particular in the tetrahedral 
models the two-quantum levels split in the same 
fashion as for a central field, viz., into one 
isolated term (the 2s wave function) and one 
triply degenerate term (three arbitrary orthogonal 
linear combinations of the usual 2p wave 
functions). It is to be particularly noted that 
Mulliken’s tetrahedral model does not demand or 
even allow “hybridization,” i.e., intermingling of 
the 2s and 2 wave functions in linear combi- 
nations. The method of calculation with the 
Mulliken model is essentially a generalization of 
one devised by Hund for simple cases, and is at 
least superficially very different from the Heitler- 
London method, as we shall see. The new contri- 
bution of Mulliken’s recent paper, as compared 
to Hund’s previous work, is the qualitative 
consideration of the arrangement of energy levels 
under molecular fields of given symmetry types. 

Both the Slater-Pauling and Hund-Mulliken 
schemes predict a right-angled model for H2O 
if one neglects (a) the forces between the two H 

*F. Hund, Zeits. f. Physik 73, 1, 565; 74, 429 (1931-2). 


7R.S. Mulliken, Phys. Rev. 40, 55; 41, 49 (1932). 
§H. Bethe, Ann. d. Physik 3, 133 (1929). 
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atoms, (b) the finite separation of the 'S, *P, 'D 
terms arising from the configuration p* of O, and 
(c) the bonding power of the two 2s electrons of 
O. This right-angled feature has already been 
demonstrated common to both methods in the 
literature, but we may here call attention to the 
fact that they give nearly similar energies of 
formation of HO when the constants are 
determined by an assumed value of the energy of 
OH. Namely both schemes are readily seen to 
yield the result that the energy of the reaction 
2H+0=H.0O is double that of H+O=OH, if we 
neglect items (a), (b), and (c) above, also (d) 
interaction terms between 1s of H and 27 of O, 
which are surely subordinate to those between 
1s and 20, and (e) allowance in the Heitler- 
London-Slater-Pauling method for the fact that 
the electron pairing is imperfect, i.e., the valences 
not 100 percent ‘‘pure’’ in Weyl’s sense. Slater 
mentions that the heat of formation of OH from 
O and H is probably about 5 volts, and hence 
both schemes agree in predicting that the heat of 
H.O from O and 2H is ia the neighborhood of 10 
volts, in accord with experiment.’ 


CRITICAL COMPARISON OF MULLIKEN AND 
SLATER-PAULING WAVE FUNCTIONS: PROS 
AND CONS OF HYBRIDIZATION 


Our next point will be a mathematical calcu- 
lation exhibiting explicitly why the s and p wave 
functions do not hybridize in the Mulliken model 
of CH,. We shall start with a system of repre- 
sentation in which the sp* wave functions are 
used, and show that proper solution of the 
secular problem introduces a canonical transfor- 
mation which destroys the hybridization. Fol- 
lowing Hund, we shall assume that an idealized 
localized valence bond between an sp* wave 
function of C and a 1s wave function of H can be 
represented by using a wave function of the form 
ay(1s; H)+(1—a?)!¥(sp?; C) where the axis of 
the sp* function is, of course, to be taken co- 
incident with the C—H axis. When we use 4 
notation such as ¥(1s; H;) we mean a 1s wave 


9 A. S. Coolidge has recently made a calculation of H:0 
by the Slater-Pauling method without the use of as many 
approximations as (a)—(e); (Phys. Rev. 42, 189 (1932)). 
Unfortunately the accord with experiment is poor, sug- 
gesting that the close agreement achieved with the less 
refined calculations is perhaps somewhat accidental. 
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function for an H attracting center. The sub- 
script 7 (t=a, b, c, d) is inserted only when it is 
necessary to distinguish between the four differ- 
ent H nuclei involved in CH,. Similarly ¥(2p; C) 
denotes a 2p wave function for a C attracting 
center, ¥(sp*; C) an sp* one, etc. When wecome to 


I=y¥(1s;H.), I=y(is; Hy), M1l= (1s; H-) 


The secular determinant is of the form 





a—W b b b 

b a—W b b 

b b a—W b 

b b b a—W 
c—dW e-—fW e-fW e—W 
e—fW c—dW e-—fW e—-fW 
e—fW e-—fW c-—dW e-fW 
e—fW e-—fW e-—fW c—dW 


The rows and columns are indexed with respect 
tol’, ---, IV’, I, ---, IV. The constants c, e are 
Hund resonance terms between the C and H 
atoms, and g a similar term between the H atoms 
themselves. By a resonance term or integral we 
mean a matrix element /¥/Hydv in which y’, y 
are single electron wave functions for different 
atoms and in which H is the Hamiltonian 
operator. It is not to be confused with ‘“‘exchange 
integrals,” which have a meaning only in two-or- 
more electron systems and in which electrons are 
permuted rather than simply one electron trans- 
ferred from one atom to another. The term } 
exists because the hybridization prevents / being 
a good quantum number. The energy W appears 
off the diagonal because the wave functions of 
different attracting centers are not mutually 
orthogonal. On the other hand the C wave 
functions are orthogonal to each other, and hence 
W does not appear off the diagonal in the upper 
left quarter of (2). We shall not attempt at 
present to express the various constants in terms 
of the familiar integrals of the two-center 
problem, but substantially this will be done in 
Part II in connection with more general models 
(cf. especially note 7 of Part II). The occurrence 








CH, we can no longer hope to express the 
molecular wave functions as linear combinations 
of two wave functions, as in a diatomic problem, 
but instead must use linear combinations of 


the following eight wave functions, denoted by 
I-IV’ 


IV=y(1s; Ha); 


I’, II’, III’, 1V’ =sp* wave functions of C directed along 


(1) 
C—H,, C—H,, C—H., C—Haz respectively. 
c-dW e—fW e-fW e-—fW 
e—{W c—dW e-fW e-fW 
e—fW e-—fW c—dW e-—fW 
e—f{W e—-—fW e-fW c—dW 
k—W g-hW g-hW g-hWw @) 
g—-hW k-W g-hW g-hW 
g—-hW g-hW k—-W- g-hw 
g-hW g-hW g-hW k-—-W 





of the same letter in several places in (2) is, of 
course, a consequence of the tetrahedral sym- 
metry and complete parity of the four sp* wave 
functions except for direction. Clearly the Hund 
resonance term has one value c for pairs of wave 
functions such as (I, I’), (II, II’)--- (IV, IV’) and 
another common value, say e, for (I, II’), 
(II, I’), (1, III’), etc., while g denotes the value 
for H—H pairs such as (I, II), ete. 

The fact that the nonvanishing off-diagonal 
elements of (2) are not confined to the places 
(I; 1’), (II; Il’), ---, (IV; IV’), of course, shows 
that localization is impossible, i.e., that twofold 
linear combinations of the form al+(1—a?)!I’, 
etc., are inadequate wave functions. Great simpli- 
fication can, however, be achieved by making a 
canonical transformation corresponding to the 
introduction of the eight new wave functions!” 


10 In unpublished work, Mulliken has also independently 
found that the functions ¥, 2, ¥2, v4 defined in (3) are 
linear combinations of the H wave functions appropriate to 
the tetrahedral symmetry. Mulliken’s method of obtaining 
the combinations was a study of general symmetry prop- 
erties rather than of the secular determinant. He has also 
thus deduced independently the trigonal wave functions 
to be given in Eq. (29) of Part II. 
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Wi=$3(I+II4+II4+IV], y2=3{1+lN—-WI—-IV], 


¥s=301-HU+4+II-IV], ws=3(1-N-II+Iv] 


(3) 


vi =3( +I +1 +1V'), ye = 31’ + 1’ - Ul —- IV’), 
vs =3(1/-I+11l—-IV’], 4’ =3(l’- 11-11’ +IV’] 


in place of I---IV’. One verifies that the secular determinant then factors into four second order 
determinants, of which three are identical. The determinantal equations occurring once and thrice are 
given respectively by the upper and lower choice of sign in the following Eq. (4): 


a+(b+2b)—W 


The pairs of wave functions which are associated 
with the four second order determinants are 
respectively (Y1, ¥1'), (Wa, Wo’), (Wa, Hs’), (War Wa’). 
All but the first of these pairs are correlated with 
the lower choice of sign in (4). The transfor- 
mation (3) which simplifies (2) to (4) completely 
eliminates the sp* hybridization, since y:' does 
not combine with y2’, ¥3’, ys’, and since further 
y,’ is merely an ordinary 2s wave function. The 
latter statement follows from the fact that the sp* 
functions are of the form? + $y¥(2s)+f(r)(r-1,), 
where 1; is a unit vector directed along C—H,i, 
whence 1,+1,+1.+14=0. Also yo’, 3’, Ya’ are 
ordinary 2p wave functions. 

One thing must be particularly noted. Ac- 
cording to the ideas of Hund and Mulliken we 
must fill twice each of the lower roots of the four 
second order equations :—twice because the Pauli 
principle allows two electrons to a given orbital 
state. The total wave function for the eight 
electron CH, system is (forgetting, as we do 
throughout, the K shell of C) 


(1), (2) (3) b2(4)H3(5)b3(6)b4(7)H4(8) (5) 
with 
&;=6.yi+(1-—B2)!y7 (i=1, as 4) 
B2=B3=Bs. (6) 


Here #,(j) means that the coordinates involved 
as arguments of #, relate to electron j. The 
values of 8;, B2 are determined by finding the 
transformation matrices connected with the two 
secular problems (4). 

Accumulation of electrons on one nucleus. If we 
express (5) in terms of the original functions 
I, ---, IV’ by means of (3) we see that (5) will 
involve terms of the form I(1)I(2)---1(8). The 


c+(e+2e) -—dW—(f+2/)W) 
c+(e+2e)-—dW—-(f42f)W k+(g+2¢)-—W-(h42h)W] | 





(4) 


physical interpretation of such terms is that 
there is a finite probability of all eight electrons 
being piled momentarily on one hydrogen atom, 
since I(j7) means a hydrogen wave function in 
which electron j is attached to nucleus a. There 
are also other terms involving accumulation of 
from 2 to 7 electrons on the same H atom. 
Analogous results are found in Hund’s calcu- 
lations for H2 and H2O, where respectively 2 and 
4 electrons can accumulate on an H atom. 
Fortunately, all these remarks apply only as 
long as we take the wave functions to be simple 
products such as (5) of ‘‘one-electron functions.” 
Actually, the seriousness of the accumulation 
effects is considerably diminished from the 
foregoing when one includes the spin and “‘anti- 
symmetricizes” so as to satisfy the Pauli princi- 
ple. One must, of course, take the antisym- 
metrical linear combinations, inclusive of spin 
functions, since the wave function (5) is not 
unique because of the possibility of permutation 
of electrons. As long as we have agreed to include 
electron repulsions only insofar as they can be 
represented by the self-consistent field, the 
process of taking such combinations will have no 
effect upon the energy provided we neglect terms 
resulting from the want of orthogonality of the 
wave functions I, I’, ---, IV’ etc., i.e., provided 
we neglect the elements dW, fW, hW in (2). 
When linear combinations are taken in ac- 
cordance with the Pauli principle, it is found that 
not more than two electrons can accumulate on a 
given nucleus. In other words, the final wave 
function will involve terms such as I(1)1(2)I1(3) 
-++, but none such as I(1)1I(2)1(3)---. For this 
observation the writer is indebted to Professor 
Slater. Similarly in Hund’s model of HO, only 
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two electrons can accumulate on one H nucleus 
when the exclusion principle is taken into 
account. 

To prove that the exclusion principle reduces 
accumulation as stated above, one utilizes the 
well-known result of Slater’s paper on complex 
spectra," that the complete wave function, in- 
clusive of spin, can be written as a sum of 
determinants. The elements of a determinant 
consist of the product of a spin factor and an 
orbital wave function. Elements of the same 
column have the same coordinate and spin 
arguments, i.e., relate to the same electron, but 
have different spin or orbital functions, while 
elements of the same row all involve the same 
functions but with different arguments. In our 
case the orbital factors of each element may 
without loss of generality be considered to be one 
of the eight expressions I---IV’, since provision 
for linear combinations is made when it is stated 
that the complete wave function is the sum of 
determinants. Because the spin factor can only 
assume two different values 6(—4; ms), (3; ms), 
any given one of the expressions I---IV’ cannot 
appear in more than two rows of any determinant 
without the latter vanishing. From this it follows 
that the expansion of the complete wave function 
cannot involve terms in which any one of these 
expressions appears as a factor more than twice, 
Q.E.D. 

Another unpleasant feature is that much of the 
expansion of (5) represents terms of the struc- 
tures C4*(H,)*-, C3+(H,)*-, ---, C4-(H4)**, even 
though in the mean the molecule is not ionic. For 
instance, if 6,2=6.2?=4 and if we neglect correc- 
tions for non-orthogonality, the C and collec- 
tively the four H atoms share the 8 electrons 
equally on the average but one calculates that 
then the molecule spends the following per- 
centages of the time in various structures: 
C%(H,)°, 27 percent; C+(H4)~, 22; C-(H,4)t, 22; 
C+(H,)?, 11; C*(H,)**, 11; C#(H,)*, 3; 
C*(H,)+, 3; C4+(Hy)*, 3; C*(H,)**, 3. (The 
process of antisymmetricizing does not affect the 
percentages of distribution between C and 4H, 
but does show that in C‘+(H,)* the 8 electrons 
are always distributed evenly between the 4 
hydrogen nuclei.) This is entirely too much in- 
stantaneous polarity. The excessive instantane- 


Titties 


"J. C. Slater, Phys. Rev. 34, 1293 (1929). 


ous polarity Cemanded by the Hund-Mulliken 
method is glaringly revealed in the case of He, 
where this method predicts the molecule to be 
H+H~- and H-H¢* half the time, whereas calcula- 
tion shows that it probably is so constituted only 
1/64 of the time.” If the transfer effects are as 
frequent in CH as in Hp, and if those in CH, can 
be roughly estimated by fourfold superposition 
of those in CH, then this value 1/64 suggests 
that methane should be C+(H,)~ or C~(H,)* 
each only 3 percent of the time and be in doubly 
or more ionic states (all kinds combined) less 
than 1 percent. The excessive accumulation 
effects can be eliminated by extending Mulliken’s 
method by adding to (5) other solutions in which 
one or more electrons settle in the higher’ roots 
of the determinants (4). Essentially because the 
presence of nodes increases materially the nor- 
malization constants for the higher states, the 
amount of such states which must be added to 
(5) is surprisingly small, at least in He, where 
Mulliken notes” that the accumulation effect is 
properly reduced by adding’ y,(1)¥.(2) to 
¥,(1)¥,(2) in the ratio of only 1 to 8. This 
suggests that it probably would be quite a good 
approximation if we could make a calculation in 
which (5) is augmented by linear combination 
with states in which only one or two electrons 
settle in the higher roots of (4). Such a calcula- 
tion, however, presents serious mathematical 
difficulties, and without it the Hund-Mulliken 
method must be regarded more as a qualitative 
rather than quantitative procedure. 

The preceding result that two electrons can 
occasionally accumulate on one atom is, of 
course, in sharp contrast to the high degree of 
localization found in the Slater-Pauling model, 
and related Heitler-London methods of calcula- 
tion which assign one and only one electron to an 
H atom. Naturally, accurate inclusion of the 
electron repulsion prevents accumulation being 
realized to the extent predicted by the preceding 


12 See discussion by R. S. Mulliken, Phys. Rev. 41, 69-71 
(1932) based on a calculation by Slater, ibid. 35, 514-515 
(1930). 

138 Do not confuse our references to the higher and lower 
roots of an individual determinant of (4) with the upper 
and lower choices of sign which yield the two determinants. 

14 Here, as usual, the notation yy, ¥, means wave func- 
tions respectively even and odd with respect to reflection in 
the midpoint of the internuclear axis. The corresponding 
notation in reference 12 is go, ¢1. 
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calculation. We have agreed to neglect this 
repulsion except via the self-consistent field, and 
under such neglect the question of accumulation 
is irrevelant for the time being, but any wave 
functions giving pronounced accumulation will 
be rather severely perturbed, and hence cannot 
be regarded a good approximation, if the theory 
is refined by introducing accurate consideration 
of the electron repulsions as a perturbation. 
Strict localization of one electron to an H atom 
could have been achieved if we had not taken the 
linear combinations (3) but had instead kept the 
hybridized sp* wave functions and had used 
only twofold linear combinations of the form 
al+(1—a?*)!I’ etc., which, as already stated, 
represent purely localized bonds. Such a mode 
of securing localization would involve the error of 
neglecting the elements J, e, f, g, h in (2) since 
these elements cannot be brought onto the 
diagonal by taking such twofold (rather than 
eightfold) linear combinations. This error we 
shall call for brevity the neglect of individual 
symmetry type inasmuch as the wave functions 
representing sp* bonds are thus inadequate to 
diagonalize (2) or in other words lack the 
symmetry properties characteristic of irreducible 
representations of the tetrahedral group. The 
word “individual’”’ is inserted because the sym- 
metry considerations relate to wave functions 
of an isolated electron. When electron repulsions 
are accurately included, the dynamical problem 
ceases to be that of a one-electron system, and so 
symmetry properties for individual electrons are 
no longer preserved. Nevertheless, they will not 
be seriously impaired if the Hartree self-consist- 
ent field has been a good approximation. 

The seriousness of neglecting symmetry type is 
gauged by the magnitude of B, e, f in (2), or more 
succinctly, by the separation of the lower roots” 
of the two determinants (4) since if these roots 
coincided their wave functions could be combined 
by a process more or less the reverse of (3) into 
those representing localized sp* bonds. If the 
carbon levels are lower than the hydrogen ones, 
the two lower roots of (4) are W=a+30, 
W=a-—b when c=d=e=f=0, as is the case 
when all interaction between the C and H atoms 
is neglected. Hence 4b has the significance of 
being the carbon s-p separation, and neglect of 
individual symmetry type will be a poor approxi- 
mation if this separation is greater than or 
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comparable with the C—H _ bonding energy. 
Actually the s-p interval in carbon amounts to a 
few volts, and hence is not negligible compared to 
the energy of a C—H bond (4 volts). (The value 
of the s-p interval in C depends on which 
spectral terms are selected from the configura- 
tions sp* and s*p*; the lowest state *P of s*p* is 
1.6 volts below °S of sp’, according to rather 
indirect evidence from band spectra," and is 
definitely known to be from 8 to 10 volts below 
most’ other terms ('P,°P,'D,*D) belonging to sp’. 
Since a one-electron model is inadequate to take 
account of different term values for the same 
configuration, it is perhaps fairest for our 
purposes to compare the mean of all terms arising 
from s*f? with the mean of all terms from s?*. 
When allowance is made for the fact that the two 
other terms 'D, |S of s*p? are 1 and 2 volts higher 
than *P, the difference between the two means is 
found to amount to 6 or 7 volts.) Even if the 
s—p separation of carbon were zero, the neglect 
of individual symmetry type would not be a good 
approximation as the two determinants (4) are 
not identical even when 6=0. The situation is 
thus quite different from that which Hund‘ 
encounters in a strictly right-angled model of 
H.O, where only » wave functions are used for 
the O valences, and where he shows that the two 
lower roots of his two determinants exactly 
coincide, so that the error due to neglect of 
individual symmetry type in securing localized 
bonds is nil. 

To summarize, the Slater-Pauling localized 
bonds avoid excessive electronic accumulation on 
the same atom, but give rather too great de- 
partures from individual symmetry type. These 
properties were, in fact, already known in a very 
general way, but it seems instructive to us to 
have deduced them more concretely by means of 
the secular determinant (2). 

In Part II we shall give calculations which aim 
to show that the tetrahedral model of CH; is 
firmer than other models in both the Hund- 
Mulliken and Slater-Pauling methods, even 
without invoking H—H repulsions. 

The writer is indebted to Professor R. S. 
Mulliken for helpful discussion and correspond- 
ence. 

16 W. Heitler and G. Herzberg, Zeits. f. Physik 53, 5? 


(1929). 
16 The location of the remaining term sp* *S is unknow®. 
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Effects of Pressure on the Visible Band Spectra of Nitrogen 


HERBERT S. HARNED AND EUGENE R. BROWNSCOMBE,* Department of Chemistry, Yale University 
(Received January 7, 1933) 


intensity up to a pressure of 0.2 mm. At still higher pres- 
sures the second positive bands were found to fade until 
at pressures of 4 mm they were just noticeable. (3) In the 
case of the first negative bands, it is found that increase of 
pressure increases the number of N2* ions in the higher 
rotational states. 








e 
h 
bs 
$ 
i 
's (1) The visible band spectrum of nitrogen produced by 
. an electrodeless discharge has been studied at pressures 
. ranging from 0.002 to 4 mm. (2) At the lowest pressures the 
re spectrum was almost entirely composed of the first negative 
ce bands due to N2*. The first positive bands appeared very 
1e faintly. As the pressure of the discharge was increased, the 
ir first and second positive band systems greatly increased in 
Ig 
F. 
” EXPERIMENTAL 
er 
" HE spectra were excited by means of an 
be electrodeless discharge. A high-frequency 
ot oscillation was maintained in a copper helix in 
ol the center of which was the discharge tube. The 
we oscillator circuit was similar to that recently 
‘. described by Hunt and Schumb.! The exciting 
id’ coil was made of 1/8” copper tubing wound in a 
of helix 4” in diameter, the turns being 1/2” apart. 
for A 50-watt oscillator (De Forest Type 503A) was 
_ used. When adjustments were made, the plate 
tly current was 60 m.a. when there was no discharge 
al in the tube. This was convenient because when 
zed the glow in the discharge tube was highest the 

plate current rose to 100 m.a. 
zed The wave-length of the oscillator was 16 
)on meters (19X10° cycles). There is every reason 
de- to believe that the excitation was electrostatic 
ese rather than electromagnetic.” The oscillator pro- 
ery duced a continuous wave, -so that the rate of 
' oat change of magnetic flux was small. Uniform glow 
s of was obtained throughout the tube, and no 

indication of ring discharge was noticeable. The 
aim pressure range over which the discharge could be 
i, is obtained was very considerable. From 20 mm 
ind- ae ; ae ai 

rhe material contained in this contribution constituted 

es part of a dissertation presented by Eugene R. Browns- 

combe to the Graduate School of Yale University in partial 
mm’ fulfillment of the requirements for the degree of Doctor of 
ond- Philosophy, June, 1932. 

‘Hunt and Schumb, J. Am. Chem. Soc. 52, 3152 (1930). 
* MacKinnon, Phil. Mag. 8, 605 (1909); Brasfield, Phys. 

3, 52 Rev. 35, 1073 (1930); Lynch and Hilberry, ibid. 37, 1091 


(1931); Knipp and Knipp, ibid. 38, 948 (1931). 
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down to the lowest pressures obtainable with 
diffusion pumps and a liquid air trap, certainly 
less than 10-* mm, light glows were obtained. At 
the lowest pressures, the gray glow formed when 
the discharge was turned on, disappeared in a 
few seconds and gave way to brilliant variegated 
colors on the inner surface of the tube. 

The discharge tube was of quartz and was 12” 
long and 3/4” in diameter. A quartz optical 
window was fused on an end. A quartz tube is 
very desirable even when working with visible 
spectra since it can be heated in vacuo to high 
temperatures to remove adsorbed impurities. It 
was connected to the rest of the apparatus by a 
quartz to Pyrex seal. A liquid-air trap was 
inserted between the tube and the rest of the 
apparatus which served to keep mercury from 
the diffusion pumps and other impurities from 
the tube. The gas was admitted to the tube 
through this trap. 

A large Universal Schmidt and Haensch 
Spectrograph was employed. For visible regions 
the optical system contained two Rutherford 
prisms, dispersion (C—F 6° 52’). The slit 
mechanism of the spectroscope was placed close 
to the quartz window of the tube without an 
intervening lens. 

Tank nitrogen was used in most of the 
observations. This was found justifiable since the 
spectrum obtained was identical with that ob- 
tained by nitrogen prepared by heating potas- 
sium nitrite. Spectra were obtained at pressures 
of 0.002, 0.007, 0.013, 0.025, 0.05, 0.13, 0.2, 0.4, 
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Fic. 1. Effect of pressure on the band spectrum of nitrogen. 


0.8, 2.2 and 4 mm. Five-minute exposures were 
used. 


RESULTS AND DISCUSSION 


The effects of pressure on the entire visible 
band spectrum are illustrated by Fig. 1, in which 
photographs obtained at the pressures indicated 
at the right are shown.* At the lowest pressure, 
every line observable at wave-lengths less than 
5500A is due to the spectrum of the first negative 
group due to Net with the exception of the two 
lines of atomic nitrogen indicated at 5463A and 
4358A. The vibration-band heads of this system 
of bands are indicated and designated by the 
vibrational quantum numbers given at the 
bottom of the illustration. At the red end, the 
first positive group of bands due to neutral 
nitrogen appear faintly. As the pressure is 
increased, the first positive group increases 
greatly in intensity, and there also appears at 
pressures above 0.01 mm another band system 
due to neutral nitrogen and designated the second 
positive group. The vibration-band heads of the 
second positive group are marked and designated 
by their vibrational quantum numbers at the 
top of the illustration. Upon closer observation 


’For a general discussion of these band systems, see 
Weizel, Bandenspektren, pp. 346-361, Handbuch der ex- 
perimental Physik, Erganzungswerk, Band 1, Akad. Ver- 
lagsgesellschaft, Leipzig, 1931. 


it is interesting to note that the intensity of the 
second positive group increases up to a pressure 
of 0.2 mm and then decreases with further in- 
crease in pressure until at a pressure of 4 mm it is 
just noticeable. This particular point is more 
clearly indicated in Fig. 2 in which an enlarge- 
ment of that part of the spectrum in Fig. 1 which 
includes the (0, 1) and (1, 2) bands of the first 
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Fic. 2. The (0, 1) and (1, 2) first negative bands and the 
(1, 5) band of the second positive. 
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negative group and the (1, 5) band of the second 
positive group is shown. 

In discharges of this kind, it is to be expected 
that the negative bands or atomic lines persist at 
all pressures, since the presence of both positive 
ions and electrons are necessary in order for the 
discharge to persist. In similar experiments with 
oxygen, the second negative bands due to O2* 
were found to be intense over a large pressure 
range. In the case of hydrogen, it is well known 
that no spectrum due to H2* appears, but in this 
case the atomic lines are very intense at all 
pressures as are those due to the band spectrum 
of He. Thus, the presence of bands due to N2* is 
to be expected. Further, an activated neutral 
molecule can result from a three-body collision 
between a positive ion, a molecule and an 
electron, so that at the higher pressures we may 
expect the appearance of the positive groups of 
bands. The fact that the second positive group 
decreases in intensity at the higher pressures 
shows that the energy exchanges bear a compli- 
cated relation to the number of collisions. Under 
the conditions of these experiments, the energy 
contributed to the gas in producing the discharge 
was nearly the same at all pressures. Similar 
results were obtained in the region of the (0, 0) 
band of N+ at 3195A. 

The structure of the first negative bands has 
been studied by Fassbender* and Herzberg‘ and 
is of the *2—*Y type. The rotation structure 
possesses a P and an R branch but no Q branch. 
The well-known phenomenon of alternating in- 
tensities in these lines appears very clearly. 
Measurements of the lines of the N2+ bands from 
the plates obtained at 0.007 mm and 4 mm were 
made and were found to agree within +0.1 of a 
wave number. These measurements with vacuum 
correction also agreed within the same limits with 
the values given by Fassbender. The equations 
for the two branches, 


P(K) = 23391.01 —3.99K +0.19K? 
and 


P(R) = 23391.01+3.99(K +1)+0.19(K+1)?, 


were found to represent the measurements to 
within +0.05 of a wave number. From this, the 


‘Fassbender, Zeits. f. Physik 30, 73 (1924); Herzberg, 
Ann. d. Physik 86, 189 (1928). 
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moments of inertia of the lower and higher 
vibrational states were found to be 


I’ =(13.23+0.05) X 10-” g-cm? 
and 
I’ =(14.56+0.05) X 10-” g-cm?. 


13.35 and 14.41 have been given for I’ and I”’ of 
the (0, 0) band.’ The origin of the band and some 
of the rotational quantum numbers of the P and 
R branches are indicated at the bottom of Fig. 2. 
The structure of the second positive band is a 
3]I — JIT type.* No attempt was made to measure 
the lines of this band. One further fact was 
observed from microphotometer traces obtained 
at the higher and lower pressures. At higher 
pressures, the intensity of the lines at higher 
quantum numbers was relatively greater than 
their intensities at lower pressure. This is illus- 
trated clearly in Fig. 3 where the intensities, 
roughly determined from the microphotometer 
traces, are plotted against the rotational quan- 
tum numbers, in this case the even numbers, at 
pressures of 0.007 and 2.2 mm. The odd-num- 
bered lines show an exactly similar behavior. 
Since between the quantum numbers 4 and 10 
the N,* lines are augmented in intensity by the 
lines of the second positive band, it is roughly 
estimated that the distribution of intensities is 
more truly represented by the solid line than by 
the dotted one given by the measurement. 
There seems to be little doubt that pressure 
increases the number of No»* ions in the higher 
rotational quantum states. 





L0.4.0.5.4.3 0.8.2 2 
0246 10 14 18 22 
K 


Fic. 3. Intensities of the lines of (0, 1) first negative band 
of nitrogen. 
5 International Critical Tables 5, 415 (1929). McGraw- 
Hill Book Company. 
6 Lindau, Zeits. f. Physik 25, 247 (1924); 26, 343 (1924); 
30, 187 (1924). 
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The Exchange of Energy Between Polyatomic Molecules and a Metallic Surface 


FRED R. WHALEY,* Department of Chemistry, Johns Hopkins University 


Experimental evidence is given to show that fairly 
complicated molecules, on striking a hot surface, do not 
absorb energy in their vibrational degrees of freedom to any 


(Received January 13, 1933) 


appreciable extent. In regard to the energy interchange the 
molecules behave like hypothetical gases of six degrees of 


freedom, translational and rotational. 





OST of the work reported on energy ex- 

changes between gases and solids deals 
with light monatomic or diatomic molecules. 
However, Rice and Byck! in attempting to 
decompose acetone and dimethyl mercury by 
impinging it on a platinum surface at 1600°C 
concluded that the “transfer of energy from a 
heated target to the internal degrees of freedom 
of an organic molecule is not 100 percent 
efficient.”” The purpose of this work is to show 
experimentally that, when fairly complicated 
molecules strike a hot surface, the efficiency of 
energy transfer from the solid surface to the 
vibrational degrees of freedom of the impinging 
molecule is very low and, within the limits of 
experimental error, possibly equal to zero. 

The apparatus consisted essentially of a 
water-cooled, cylindrical chamber about 10 cm 
long and 4 cm in diameter, in which was sus- 
pended a platinum wire. With the aid of tungsten 
leads sealed through the glass, the wire could be 
heated electrically and its temperature calculated 
from its resistance. Allowance was made for the 
almost negligible resistance of the leads, and 
this correction was kept constant by keeping the 
leads at the temperature of solid CO» throughout 
all the readings. A bulb containing the substance 
being used and surrounded by a bath of alcohol 
and solid CO, was connected with the chamber. 
The pressure in the chamber could be measured 
by means of a McLeod gauge, so constructed 
that when pressure measurements were taken 
the pressure of the vapor in the capillary was 
only a fraction of the vapor pressure at room 


* Du Pont Fellow. 
1F, O. Rice and H. Byck, Proc. Roy. Soc. A132, 50 
(1931). 
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temperature of the substance used. The chamber 
could be evacuated either rapidly through a wide 
tube or slowly through a capillary. 

A calibration curve was obtained by evacu- 
ating the chamber to about 10-* mm of Hg and 
plotting the energy input in the wire against the 
temperature of the wire. Due to varying con- 
ditions of the surface of the chamber, probably 
because of platinum and carbon deposited there, 
the calibration curves varied slightly from day to 
day. Curves obtained before and after a run 
agreed, however, so separate calibration curves 
were made with each run. 

The run itself was made by keeping the sub- 
stance in the chamber at a known pressure and 
again plotting energy input against temperature. 
At a given temperature the energy given up to 
the colliding molecules could at once be read 
from the graph. This was compared with the 
energy exchange calculated, assuming the gas 
to be a hypothetical gas of only six degrees of 
freedom, translational and rotational. 

In this calculation the following equation 
derived by Langmuir? was used: 


m=p(M/2xRT)}, 
where m=grams striking surface per second per 
cm?; p=pressure; M=molecular weight; T =ab- 
solute temperature; R=gas constant. The mols 
per second hitting wire=mA/M where A is the 
area of the wire. The energy absorbed per second 
in watts is: 


(4.185mA(T2—-T))Cv]/M, 


where (T2—T;) is temperature difference be- 
tween wire and walls and C, is the specific heat. 


21. Langmuir, Phys. Rev. 2, 329 (1913). 
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EXCHANGE OF ENERGY 

Runs were made with chloroform, acetone, 
methyl alcohol, carbon tetrachloride and ethyl 
acetate. The substance was kept at constant 
pressure during a run by the constant-temper- 
ature bath and the slow pumping through the 
capillary. The pressure should be sufficiently 
low that the temperature distribution of the 
molecules colliding with the wire is equal to the 
distribution of molecules in temperature equi- 
librium with the walls. This means that the mean 
free path must be at least as great as the dia- 
meter of the wire, and preferably two or three 
times as great. 

In the runs with methyl alcohol, carbon 
tetrachloride and ethyl acetate the pressure was 
sufficiently low, the highest being 0.047 mm for 
ethyl acetate. The mean free path was 0.034 cm. 
The diameter of the wire was 0.025 cm. 

Let us designate the molecules about to strike 
the wire as A molecules. These had their last 
collision with other molecules, which we shall 
call B molecules, when they were at an average 
distance from the wire equal to A, the mean free 
path. The A molecules will have a temperature 
distribution equal to that of the main body of the 
gas, unless an appreciable fraction of the B 
molecules just came from the wire, hence having 
a temperature increment. That fraction is 
approximately the angle in radians subtended 
by the wire from a point at distance A, divided 
by 27. 

In the ethyl-acetate run this fraction is about 
0.12, but the A molecules, colliding with B 
molecules coming from the wire, can at best 
receive a temperature increment of half the 
increment received by the B molecules. Since the 
efficiency of energy transfer on collision may not 
be perfect, the error due to this preheating of the 
A molecules is probably even less than the 
maximum error, which is 0.12/2=0.06 or 6 
percent. This is within the limits of experimental 
error. For methyl alcohol and carbon tetra- 
chloride the error is even less. Figs. 3, 4 and 5 
show the marked agreement between the ob- 
served energy changes and those calculated for a 
hypothetical gas of six degrees of freedom. 

In the chloroform and acetone runs this error 
is about 13 percent, which is too high even 
allowing for inefficient energy exchanges between 
molecules A and B. This would tend to make the 
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Fic. 1. Chloroform. 


calculated energy exchange too high. Figs. 1 and 
2 show a calculated value higher than the 
measured value. 

In each graph, 1 is the calibration curve, 2 is 
the observed energy curve and 3 is the energy 
curve calculated on the basis of perfect efficiency 
in energy transfer to translational and rotational 
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degrees of freedom and no transfer of energy to 
the vibrational degrees of freedom. Curve 4 is 
calculated by assuming perfect efficiency of 
energy transfer to all degrees of freedom. 
(Specific heat values were calculated by Dr. 
D. H. Andrews and his students in this labora- 
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tory from Raman data.) Experimental points on 
curve 2 are shown. Table I summarizes the con- 
ditions under which the experiments were carried 


out for each of the substances studied. 


TABLE I. 








Ratio of 
A to 
wire 

diam. 


0.56 
0.52 
6.10 
2.64 
1.36 


A at 0°C 
and one 
atm. 
(cm) 


265 X10-* 
258 x10-* 
517 X10-% 
326 X10-% 
212 x10 


Pressure 
(mm) 


0.140 
0.146 
0.0256 
0.0376 
0.047 


Substance 


CHCl; 
CH;COCHs; 
CH;0H 


CCl, 
CH;COOC;H; 











DISCUSSION 


Zener* concludes from theoretical consider- 
ations that light molecules of light atoms collid- 
ing with each other exchange rotational quantum 
numbers freely, but that the transfer of vibra- 
tional energy may be difficult even in cases of 
exact resonance. Although in the work reported 
here one of the colliding bodies is platinum and 
not a “light molecule,” the results seem to 
conform with Zener’s conclusion. 


*C. Zener, Phys. Rev. 37, 556 (1931). 
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Herzfeld and Géppert Mayer*t showed from 
theoretical considerations that the time necessary 
for an organic molecule, adsorbed on a hot 
surface, to absorb energy in its vibrational 
degrees of freedom was longer than the collision 
time. The results reported here are in accord with 
that theory. 

Van Wijk® by optical methods shows that the 
accommodation coefficient for rotational energy 
in the case of He is very small. Knudsen* by 

‘ Herzfeld and Géppert Mayer, Zeits. f. physik. Chemie 
Bodenstein Festband, 669 (1931). 

° Van Wijk, Zeits. f. Physik 75, 584 (1932). 

6 M. Knudsen, Ann. d. Physik 6, 129 (1930). 


measuring radiometer pressures finds that the 
accommodation coefficient for the internal energy 
of Hz can be equated to the accommodation 
coefficient for the translatory energy. The work 
of both these authors with He seems to be con- 
trary to the results reported here, but the latter 
deal with heavier and more complicated mole- 
cules. 

The author wishes to express his sincere appre- 
ciation to Professors D. H. Andrews and K. F. 
Herzfeld and Dr. J. E. Mayer for their kind 
suggestions and advice during the carrying out 
of this research. 
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The dipole moments of four fluorochloromethanes are 
measured in the vapor state. The inductive effects of 
carbon-halogen bond moments are calculated and, through 
these, the moments of methyl fluoride and methy] chloride 
are used to calculate approximate values for CH2Cl2, 
CHCl;, and CFCl;, which are in satisfactory agreement 
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VOLUME 1 


with the observed values. The moments of these molecules 
are then used to calculate values for other substituted 
methanes which are in excellent agreement with the 
observed values. The value of a bond moment is treated as 
a function of its environment. 








N an earlier paper! experimental determi- 
nations of the dipole moments of di- and 
tribromo- and iodomethanes were reported and 
the results were confirmed by the measurements 
of Miiller and Sack.? In this first paper on the 
substituted methanes the moment of each mole- 
cule was treated as the resultant of the moments 
of dipoles lying in the lines joining the carbon and 
halogen nuclei. The observed moments were 
smaller than those calculated on the assumption 
that these bond moments did not affect one 
another and acted at an angle of 109.5° to one 
another, the normal carbon tetrahedral valence 
angle. Kinetic theory diameters were, therefore, 
calculated for the halogen atoms and assumed to 
cause widening of the valence angles, which were 
calculated from these diameters and from the 
carbon-halogen nuclear separations estimated 
from spectroscopic and x-ray data. The distance 
between the two chlorine nuclei in methylene 
chloride was calculated to be 3.00A, while Wierl* 
later found from electron diffraction a distance of 
3.16+0.06A in methylene chloride and 3.04 
+0.04A in chloroform. In methylene bromide, 
the distance between the bromine nuclei was 
calculated to be 3.20A, while Wierl found 
3.35+0.04A in carbon tetrabromide, and, in 
methylene iodide, the distance between the 
iodine nuclei was calculated to be 3.62A, while 
Huggins and Noble,‘ by means of x-ray analysis, 


1C, P. Smyth and H. E. Rogers, J. Am. Chem. Soc. 52, 
2227 (1930). 

2H. Miiller and H. Sack, Phys. Zeits. 31, 815 (1930). 

3 R. Wierl, Ann. d. Physik [5] 8, 521 (1931). 

‘4M. L. Huggins and B. A. Noble, Am. Mineral. 16, 519 
(1931). 
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found 3.56A in crystalline iodoform and reported 
values 3.86 and 3.80+0.2A as found by Brockway 
and Pauling for iodoform in benzene solution. 
The angles between 120° and 135° obtained from 
these calculated distances gave calculated dipole 
moments in excellent agreement with the ob- 
served, an agreement which was termed some- 
what fortuitous because of the speculative nature 
of the reasoning, uncertainties in some moment 
values, and neglect of the inductive effect of 
each dipole upon the rest of the molecule. 

Smallwood and Herzfeld' showed that this 
inductive effect could lower the moment of an 
ortho-disubstituted benzene from the value calcu- 
lated as the resultant of two mutually inde- 
pendent dipoles acting at an angle of 60° to each 
other to the observed value, a not inconsiderable 
decrease. As the dipoles in the halogenated 
methanes are closer together than those in the 
ortho-disubstituted benzenes, it is evident that 
some attempt should be made to take account of 
their inductive effects. 

Wierl’s electron diffraction experiments give 
larger carbon-halogen distances than do Pauling’s 
atomic radii® and valence angles between 110° 
and 115° for the chlorinated methanes, not 
establishing the existence of any widening of the 
valence angle beyond the regular tetrahedral 
angle rounded off to 110°. In view of this, we 
shall examine the moments of chlorine-substi- 
tuted methanes with 110° valence angles assumed 
between the bonds and shall tentatively make a 


5H. M. Smallwood and K. F. Herzfeld, J. Am. Chem. 
Soc. 52, 1919 (1930). 
6 L. Pauling, Proc. Nat. Acad. Sci. 18, 293 (1932). 
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similar assumption in the case of fluorine- taining both fluorine and chlorine substituents 
substituted methanes, although calculation shows prepared by him. 
the possibility of a greater widening of the angle 


between two C—F bonds than that in the case of 
two C—Cl bonds. The dielectric constants « of the four vapors 


The kindness of Dr. Albert L. Henne of Ohio were measured with an apparatus and method 
State University has made it possible for us to described elsewhere. From each measurement of 
measure four pure samples of methanes con- the dielectric constant, the polarization was 


EXPERIMENTAL RESULTS 


TABLE I. Experimental data. 








CHFCl, CHF.Cl CFC; 
p(mm) P T(°K) p(mm) F T(°K) p(mm) eg 


426.3 51.36 303.9 971.5 55.95 299.2 586.6 29.00 
313.5 51.01 748.2 55.69 577.5 29.16 
246.0 50.85 579.4 55.50 541.3 29.23 
244.3 50.96 413.9 55.30 518.3 28.97 
227.6 50.83 273.0 55.06 424.4 28.71 

0 50.55 228.8 54.86 293.7 28.61 
157.7 54.73 261.6 28.62 
388.6 91.12 0 54.75 250.2 28.49 
335.1 50.97 0 28.06 
329.4 50.98 679.5 54.69 
283.4 50.92 674.3 54.66 499.5 28.69 
267.5 50.74 672.7 54.70 490.3 28.70 
231.2 51.00 467.4 54.46 436.6 28.72 
157.5 50.72 283.0 54.27 359.0 28.58 

0 50.55 219.6 53.90 277.3 28.27 
0 53.82 182.6 28.22 
467.5 47.63 180.7 28.20 
358.3 47.55 659.8 50.24 0 27.92 
344.1 47.45 526.6 50.05 
319.7 47.64 447.2 50.14 565.8 28.48 
247.0 47.16 283.1 49.69 545.3 28.46 
212.8 47.35 0 49.65 515.9 28.35 
184.1 46.80 515.7 28.42 

0 47.10 706.4 46.18 422.3 28.38 
586.8 46.14 369.6 28.22 
590.8 44.34 451.6 46.13 347.2 28.28 
494.4 44.18 400.7 45.91 290.5 28.12 
403.1 44.05 223.4 45.70 279.8 28.35 
382.5 44.08 0 45.80 272.9 28.19 
338.6 44.11 0 27.88 
329.1 44.08 760.0 43.23 
295.8 44.00 602.9 43.17 529.1 28.08 

0 43.70 . 591.1 42.97 510.8 28.04 
532.4 42.99 456.0 27.95 
639.6 41.95 469.4 43.28 366.1 27.82 
519.6 41.85 446.7 42.82 269.7 27.73 
474.1 41.74 0 42.80 261.6 27.72 
371.3 41.43 193.5 27.25 
357.9 41.57 721.6 40.19 189.9 27.40 

0 41.14 651.0 40.30 186.9 27.44 
612.5 40.16 182.9 27.01 
546.9 40.49 0 27.40 
455.5 40.18 
446.9 40.25 568.2 27.61 

0 40.25 541.0 27.66 
455.6 27.56 
360.6 27.48 
350.7 27.78 
264.9 27.45 
257.4 27.53 
256.8 27.27 

0 27.23 
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calculated as 
P=((e—1)/(e+2)]V, (1) 


where V is the molar volume. As the ideal gas 
law was used in calculating P, the dielectric 
constant was measured for several different 
pressures at each temperature and the calculated 
polarizations were extrapolated to zero pressure 
to obtain a value of the polarization free from 
error caused by deviations from the ideal gas 
law. It may easily be shown that, at these 
pressures, the polarization-pressure curve is 
practically a straight line with a slope pro- 
portional to the second virial coefficient of the 
gas. Extrapolation of the polarization to zero 
pressure is, therefore, a simple matter, although 
the accuracy of the polarization values is in- 
sufficient to give a satisfactory value for the 
second virial coefficient. 

Table I gives the values of the polarization for 
the different pressures at each temperature on the 
absolute scale, the extrapolated values being 
given for zero pressure. These latter values were 
used to calculate for each substance the constants 
of the Debye equation, 


P=a+0/T. (2) 


The value of the atomic polarization was 
calculated as P,4=a—Pz, the molar refraction 
for infinite wave-length, Pz, being obtained by 


TABLE II. Constants of substituted methanes. 











Per Pa b x 1018 
CHFCl, 17.0 16.1 0.9 10,260 1.29 
CHF:.Cl 15.1 11.2 3.9 12,040 1.396 
CFC; 23.8 20.8 3.0 1,256 0.45 
CF.Cl, 20.1 16.0 4.1 1,600 0.51 








TABLE III. Moments (X10"8) of substituted methanes. 








CH;F* 1.803 CCl, o> CH;Cl 1.85> 
CHFCl, 1.29 CF.Cl, 0.51 CH:Cl, = 1.57°¢ 
CHF:Cl = 1.39 CFCI; 0.45 CHCl; 1.05> 








® Measurement in this laboratory to be published shortly. 

>’Smyth, Dielectric Constant and Molecular Structure, 
Appendix I, New York, The Chemical Catalog Co. Inc., 
1931. 

eK, L. Wolf and O. Fuchs, Stereochemie (K. Freuden- 
berg), p. 244, Leipzig, Franz Deuticke, 1932. 
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extrapolation from the sums of atomic refractions 
or of electron group refractions for the sodium D 
line.? The dipole moments were calculated from 
the values of b as w=0.0127310-18(b4) and 
should be in error by no more than 0.01 X 107"8. 
The values of these constants are given in Table 
II, while Table III lists the moments together 
with those of other molecules needed for the 
discussion of the results. 


DISCUSSION OF RESULTS 


In view of the very small difference, 0.04 
10-18, between the moment of methyl fluoride 
and that of methyl chloride, the not incon- 
siderable moments of CF2Cle and CFCl; and the 
marked increase in moment caused by replace- 
ment of chlorine in CHCl; by fluorine show that 
the effects of induction and of the widening of 
valence angles by repulsion must be considered. 
The zero moments recently reported by van 
Arkel and Snoek® for the moments of CFBrs and 
CCI1;Br in benzene solution may be attributed to 
the fact that moments as small as these should be 
cannot be accurately determined by measure- 
ments on liquids. Although it was possible to 
explain the moments of CH2Cl, and CHC; as the 
resultants of moments in the C—Cl bonds, the 
angles between which were widened beyond 110° 
by mutual repulsion of the chlorines, an analo- 
gous explanation of the fluorochloromethane 
moments is possible only if the fluorine atoms 
give rise to a much smaller mutual repulsion 
than that between the chlorines. As a rough 
calculation of the internuclear distances in 
methylene fluoride carried out in the same way as 
those for the other methylene halides indicated a 
greater repulsion than in the case of the chloride, 
it would appear that the explanation is to be 
sought primarily in the inductive effects. This is 
in accord with the approach to the problem 
announced at the beginning of the paper because 
of the closeness of the valence angles of the 
chloromethanes to 110° found by Wierl and 


7 Landolt-Bérnstein, 5th edition, p. 985. Smyth, Dielec- 
tric Constant and Molecular Structure, p. 153, New York, 
The Chemical Catalog Co. Inc., 1931. 

8A. E. van Arkel and J. L. Snoek, Zeits. f. physik. 
Chemie [B] 18, 159 (1932). 

*Smyth, Dielectric Constant and Molecular Structure, 
p. 168, New York, The Chemical Catalog Co. Inc., 1931. 
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because of the considerable inductive effects 
calculated by Smallwood and Herzfeld for the 
ortho-disubstituted benzenes. 

The polarizability a of a portion of a molecule 
may, as a rough approximation, be treated as if 
the polarization occurred at one point. The 
moment induced at this point by a permanent 
moment m elsewhere in the molecule can then be 
calculated by the method of Smallwood and 
Herzfeld. If the dipole m exerting the force E is 
located at the origin of a set of plane polar 
coordinates, the horizontal and vertical com- 
ponents of E(E, and E,) at the point r, where the 
polarizability acts, are given by 


E,=m,(3 cos? 6—1)/r’+3m,(sin @ cos 0)/r 


3 
E,=3m,(sin @ cos 0)/r°+m,(3 sin? @—1)/r (3) 


where m, and m, are the components of m, which 
lies in the x, y plane. The moment induced by m 
at the point r has the components 


Miz=E,a and mj,=Ey,a. (4) 


The treatment of the effect of a system of 
charges as being that of a single dipole located at 
a point within the system becomes an extremely 
rough approximation when the effect is con- 
sidered at distances as small as those involved in 
the halogenated methanes. Although it is cer- 
tainly justifiable to treat the principal dipole of 


CH;Cl as lying in the C—Cl line, its size and its - 


location in this line are uncertain. In the case of 
the dihalogenated benzenes, Smallwood and 
Herzfeld found it necessary to locate the dipoles 
arbitrarily at the circumference of the carbon 
atoms to which the particular groups were 
attached. Consideration of the locations of the 
centers of gravity of the charges assignable tc the 
C—Cl portion of the molecule places the princi- 
pal dipole 7/8 of the distance, 1.76A,° from the 
carbon nucleus to the chlorine nucleus.!° A 
moment m located at this point would induce a 
moment in each of the three C—H bonds, which 
is calculated in terms of m by using Eqs. (3) and 
(4), the arbitrary assumption being made that a 
polarizability, a=0.66X10-*%5." acts at the 
midpoint of the C—H bond, 0.53A from the 
carbon nucleus. If the presumably small and 


‘Cf. L. Meyer, Zeits. f. physik. Chemie [B]8, 27 (1930). 
" Smyth, Phil. Mag. 50, 361 (1925). 
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unknown moment of the C—H bond is neglected, 
the observed moment, 1.85X107'’, may be re- 
garded as the sum of m and the components of the 
three induced C—H moments along the pro- 
jection of the Cl—C line. m is thus calculated as 
1.15<10-'§ and the moment induced by it in 
each C—H as 0.70X10~-'%. This value 1.15 
X10-'8 is, however, too high for the permanent 
moment of the C—Cl bond as no account has 
been taken of the inductive effect of the C—H 
moments, 0.7010-'8, back upon the C—Cl. 
This effect may be calculated and a still smaller 
value may be used for the C—Cl moment, the 
moments induced in the C — H bonds being taken 
into account in all subsequent calculations. 
However, in view of our ignorance of the moment 
of the C—H bond arising from possible inherent 
electrical unsymmetry and of other uncertainties 
which cannot be eliminated from the treatment, 
it appears rational, in cases where at least one 
hydrogen is attached to the carbon, to assign the 
entire moment of the methyl halide to the 
carbon-halogen bond. The results of such treat- 
ment will not differ greatly from those to be 
obtained by using a smaller carbon-halogen bond 
moment and adding the effects of the C—H 
bonds. 

In the methylene chloride molecule two mo- 
ments of 1.85 x 10-!8 each in the C— Cl lines at an 
angle of 110° to one another would have a 
resultant, 2.1310-18, as compared to the ob- 
served value, 1.57X10-'*. With the dipoles 
located 7/8 of the distance from the carbon to the 
chlorine nucleus, that is, 1.54A from the carbon 
nucleus, and polarizabilities, a=2.51 10-5 1 
assumed to be acting at the same points,—a 
reasonable assumption,—application of our equa- 
tions shows that the inductive effect of one 
moment of 1.85 lowers the moment in the other 
C—Cl line to 1.36. The resultant of two such 
moments happens to be identical with the ob- 
served moment of methylene chloride. However, 
the moment of 1.36 would have a smaller 
inductive effect and the lowering calculated for 
the other C—Cl moment would not be as great. 
This, in turn, would have a smaller inductive 
effect upon the first moment. Successive calcu- 
lations lead to a value 1.49 for the moment of a 
C—Cl bond when a second chlorine is attached to 
the same carbon. This would give a moment of 
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1.71 for methylene chloride if the two C—Cl 
moments acted at an angle of 110° to each other, 
but the induced moments do not lie in the C—Cl 
lines, so that the resultant C— Cl moments make 
a slightly smaller angle with each other and have 
a resultant of 1.79, which is 0.22 higher than the 
observed value for methylene chloride. 

Since one C—Cl dipole lowers the moment of 
another from 1.85 to 1.49, one might, as a very 
rough approximation, assume that two C—Cl 
dipoles would produce a lowering of 20.36 in a 
third C—Cl moment. In chloroform, which has a 
moment =3(C—Cl) cos 70°=(C—Cl), the mo- 
ment would, therefore, be 1.13 as compared to an 
observed value, 1.05. A more accurate treatment 
involving the calculation of the components along 
three reference axes leads to a value of approxi- 
mately 1.35. Obviously, if the assumed location 
of the permanent dipole is moved in toward the 
carbon nucleus, the inductive effect will be 
increased. Thus, with the dipoles and the polar- 
izabilities located 2/3 of the distance from the 
carbon to the chlorine nucleus instead of 7/8, 
still much farther from the carbon than assumed 
by Smallwood and Herzfeld, a C—Cl moment of 
1.85 reduces the other C—Cl to 1.03 instead of 
1.36 as in the previous case. Although better 
agreement between the calculated and the ob- 
served values of the moments would result from 
this moving in of the dipoles, such agreement 
would not be justified by the nature of the 
treatment. It, therefore, seems better to adhere 
to the original assumption as to the dipole 
location, an assumption which is theoretically 
plausible. Similar calculations show that because 
of its much smaller polarizability, a=0.626 
x 10-*, the C— F moment is less reduced by the 
inductive effect of an adjacent C—Cl than is a 
C—Cl moment, and because of the smaller C—F 
distance, 1.41A, the C—F moment reduces an 
adjacent C—Cl moment more than another 
C—Cl moment does. Table IV gives the lowering 
of one carbon-halogen bond moment by the 
inductive effect of another on the same carbon. 

The reason for the not inconsiderable moment 


TABLE IV. Lowering of carbon-halogen moment ( X 10'*). 
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of CFCI; is apparent in the marked differences in 
the moment lowerings in Table IV. The com- 
ponents in the F—C line of the moments induced 
by the C—F in the three C—Cl’s are so large as 
to cancel a large fraction of the components of the 
permanent moments in this line, while the C—F 
moment is less affected by the three C—Cl’s. The 
resultant difference calculated, 1.2, is much 
larger than the observed moment, 0.45. It is 
probable that the methyl halide moments as- 
signed to the carbon-halogen bonds include so 
much of the H—C bond moment, at least, part 
induced, and, perhaps, part permanent, that they 
are considerably too high for use when no 
hydrogen is attached to the carbon. The calcu- 
lation of the moments induced in the H—C’s in 
methyl! chloride gave a moment of 1.15 for the 
C—Cl and a similar calculation for methy! 
fluoride gives 0.91 for the C—F, the effects of the 
induced H—C moments upon the halogens being 
disregarded in the calculation. As this neglected 
effect raises the apparent values for the carbon- 
halogen bond moments, these may reasonably be 
rounded off to about 1.0 for C—Cl and 0.8 for 
C-—F. A more accurate calculation of these 
values is impossible because of their great de- 
pendence upon the unknown location of the point 
at which the polarizability of the H—C bond may 
be supposed to act. The use of these values in the 
calculation of the CFCl; moment gives, instead 
of 1.2, a value of 0.44 as compared to the 
observed 0.45. 

These lower carbon-halogen bond moments 
may be used equally well in calculating the 
moments of CH2Chk, CHCls, etc., because their 
inductive effects upon the H—C bonds bring the 
total resultant moments up to the same magni- 
tude as those obtained by using methyl] halide 
moments for the bond moments. It is somewhat 
disconcerting, at first, to find the use of these 
much smaller bond moments necessary in some 
cases, but consideration of the nature of a bond 
moment makes it clear that its apparent value 
cannot be constant under all conditions. The 
system of charges constituting.a carbon-halogen 
group cannot be isolated from the rest of the 
molecule, and, unless we could determine exactly 
a point in the system where its moment could be 
regarded as acting and, similarly, a location for 
its polarizability, nothing would be gained by 
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this process of isolation. A bond moment is 
merely a measure of the electrical unsymmetry of 
a certain section of a molecule and is affected by 
the environment of the section. If the immediate 
environment of the section in different molecules 
is the same, the bond moment remains constant; 
if it is different, the electrical unsymmetry of the 
section and its moment changes. For example, we 
may think of the carbon tetrachloride molecule as 
resolved into four equal C—Cl moments, which, 
because of the symmetry of the molecule, cancel 
one another to give zero resultant moment. 
Each C—Cl moment consists of the moment due 
to an inherent unsymmetry in the distribution of 
charge plus the moments induced by each of the 
three other C—Cl moments. In methyl chloride, 
the C—Cl moment contains the components in 
the C—Cl line of the moments induced by the 
H—C moments, which consist of a possible 
permanent moment plus the moment induced by 
the C—Cl moment. As only the roughest kind of 
separation of these effects can be made, it is as 
well to treat the entire moment of the molecule as 
the C—Cl moment. Attachment of more carbons 
to the methyl carbon gives molecules of some- 
what increased moment because of induction,! 
but the effect is noticeable through only two 
carbons of a chain. Little error will arise from 
treating the moment of the —CH,Cl group as a 
C—Cl moment equal to that of methyl chloride. 
In chlorobenzene, the immediate environment of 
the C—Cl is different from that in an alkyl 
chloride and the moment is somewhat different. 
If two C—Cl’s are separated from one another as 
in p- or m-dichlorobenzene, the moment of 
chlorobenzene may, with little error, be used as a 
C—Cl bond moment io calculate the resultant 
moment of the disubstituted molecule. However, 
in o-dichlorobenzene, each C — Cl bond moment is 
affected by the proximity of the other, not only 
through direct inductive effects of one C—Cl 
upon the other, but also through the effects upon 
the other adjacent parts of the molecule. The 
very approximate calculation by Smallwood and 
Herzfeld of the alteration in moment produced 
by these effects was possible because it was not 
large. From a practical point of view, the most 
useful procedure is to assign to one bond the 
moment of an entire simple molecule containing 
but one polar group, but it must be born in mind 
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that, when the immediate environment of a 
bond moment changes, the moment may change. 
The moment of CFCIl;, which has been calcu- 
lated in our previous discussion, makes an 
excellent starting point for the consideration of 
the moments in Table III. In the following 
equations, the formula of a molecule or of a 
bond enclosed in parentheses will represent its 
moment. The different bond moments are ex- 
posed to somewhat different inductive effects in 
the different molecules, but, as the total resultant 
moments of the molecules are not large, it is 
evident that the differences in inductive effects, 
which, in most cases, partially cancel one another, 
may commonly, as an approximation, be neg- 
lected. Since, however, the C—F bond in methyl 
fluoride is not exposed to the inductive action of 
other carbon-halogen bonds, one would expect 
the values calculated by Eqs. (12) and (14) to be 
lower than the moment of methyl fluoride. The 
moment of CFCI; lies in the C—F line and is 
given by 
(CFCI;) = (C—F)—3(C—Cl) cos 70° 
=(C—F)—(C-—Cl) =0.45, (5) 
since cos 70°=1/3. The moment of CF;Cl. must 
lie in the line which bisects the F—C—F angle 
and also the Cl—C—Cl angle and is given by 
(CF2Cl,) = 2(C — F) cos 55°— 2(C — Cl) cos 55° 
= 1.15[(C—F) —(C—Cl)]. (6) 
Substituting for [(C—F)—(C—Cl)] its value 
given by Eq. (5), 
(CF.Cly) = 1.15 X 0.45 =0.52, (7) 


as compared to an observed value 0.51. 

As an approximation, the moment of CHFCl, 
may be regarded as lying in the H—C line. 
Consequently, 


(CHFCI:) = (H —C) 


+(1/3)[(C—F)+2(C—Cl)]. (8) 


A similar approximation gives 
(CF:Cle) = (C—F) 
— (1/3)[(C—F) +2(C—C))]. 
Adding (8) and (9), 
(CHFCh) +(CF2Cl.) = (H— C)+(C-—F), (10) 


(9) 
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but 
(H —C)+(C—F) =(CHsF) =1.81 (obs.), 


(11) 

and 

(CHFCle) + (CF2Cls) = 1.29(obs.) 
+0.51(obs.) =1.80. (12) 


The excellence of the agreement between the 
values calculated by (11) and (12) is accidental as 
is shown by a similar calculation involving 
chloroform. 


(CHCI;) = (H —C)+(C—Cl) =1.05 (obs.). 
Adding (5) and (13), 
(CFCI;) + (CHCI;) = (H—C)+(C—F), 


or 0.45 (obs.) +1.05 (obs.) = 1.50, while (11) gives 
(H—C)+(C—F) =1.81. Subtracting (13) from 


(8) and rearranging, 
(CHFCl.) = (CHCIs) 
+(1/3)[(C—F)-—(C-—Cl)], (15) 


but (CHCI;) =1.05 and, from (5), (1/3)[(C—F) 
—(C—Cl) ]=0.15, consequently, 


(CHFCle) = 1.05+0.15 = 1.20, 


(13) 


(14) 


(16) 


as compared to the observed value, 1.29. 
As an approximation similar to that in (8), 


(CHF;Cl) = (H—C) 


+(1/3)[2(C—F)+(C—Cl)]. (17) 


AND K. B. 


McALPINE 


Subtracting (13) from (17) and rearranging, 


(CHF,Cl) = (CHCIs) 
+(2/3)[(C—F)—(C—C)], (18) 


but (CHCI;) =1.05 and, from (5), (2/3)[(C—F) 
—(C—Cl) ]=0.30, consequently, 


(CHFCl) = 1.05+0:30= 1.35, (19) 


as compared to the observed value, 1.40. 

Although the H—C moment cancels in the 
final calculations with all of these equations, it 
has been left in the initial equations to make 
clearer the assumptions. It is evident that the 
variation of moment among these different 
fluorochloromethanes agrees extremely well with 
the predictions of theory. The agreement be- 
tween the observed moments and those calcu- 
lated on the assumption of 110° carbon valence 
angles does not preclude the possibility of 
widening of the angles in, at least, some of the 
methanes. Indeed, it seems probable that there 
is more or less widening of the angles in all of the 
unsymmetrically substituted methanes, although, 
in some cases, this may be too small for detection 
by means of x-ray or electron diffraction. It 
would appear, however, that induction is much 
more important than angle widening in de- 
termining the moments of the methanes con- 
sidered in this paper. 

The writers wish to express their indebtedness 
to the National Research Council for a grant-in- 
aid, which has made this investigation possible. 
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The two-electron bond with hydrogen is found to have 
the following values of 1 x<10'8: H—As 0.10; H—P 0.36; 
H—I 0.38; H—S 0.63; H—Br 0.78; H—Cl 1.03; H—N 
1.04; H—O 1.32. These values are believed to give the 
positions of the elements on an electronegativity scale 
which is shown to be similar to that prepared by Pauling. 
The reported moments for inorganic halides as well as the 


unpublished measurements of A. L. Ferguson and M. G. 
Malone, giving ».=3.1 1078 for SbCl; and »=2.41078 
for SbBrs, give fair agreement in most cases with those 
calculated from the above scale. On the basis of calculated 
values appreciable moments are predicted for all trihalides 
of the fifth-group elements except PI; and NCI; which 
should give close to a zero moment. 





HE idea that the electric moment bears 

some relation to the ionic quality of the 
two-electron type of bond has been known for a 
number of years. It is true that the terms used 
in describing a bond have been loosely applied, 
and it is uncertain in some cases just what 
meaning has been intended. Noyes! has clearly 
described, however, the two-electron bond as 
potentially polar when it possesses a dipole 
moment. The qualitative relationship, that 
atomic bindings between like atoms do not 
possess moments while those between unlike 
atoms may do so, is generally known.’ 

The author undertook some six months ago 
the task of assigning electric moments, derived 
from the measured moments of molecules to in- 
dividual bonds. The result of this work was the 
arrangement of the elements on what Pauling® 
has since called an electronegativity scale. In 
this early attempt, the scale was built around 
carbon as a reference. The result was neither as 
consistent nor as useful as was hoped for. 

Certain qualitative regularities did appear, 
however, among the halides of the fifth-group 
elements. These will be discussed at more length 
in light of the fuller knowledge of the subject. 
The ideas herein involved were of sufficient 
interest to induce A. L. Ferguson and M. G. 


'W. A. Noyes, J. Am. Chem. Soc. 50, 2902 (1928). 

*Ruark and Urey, Atoms Molecules and Quanta, 
McGraw-Hill, 1930; Van Vleck, Electric and Magnetic 
Susceptibilities, Oxford, 1932; Glasstone, Recent Develop- 
ments in Physical Chemistry, Blakiston, 1931. 

* Linus Pauling, J. Am. Chem. Soc. 54, 3570 (1932). 


Malone to undertake experimental determina- 
tions of the electric moments of this group of 
elements. 

It was the belief of the author that any scale 
which measures the bond differences of atoms 
should relate to the scale based on electric 
moments by a constant factor. Attempts to 
prepare such a scale from independent data were 
fruitless. At this, time the paper of Pauling,’ 
giving an electronegativity scale based on 
thermal data appeared. Since the scale given by 
Pauling used hydrogen as the point of reference, 
the author revised his scale to read from the 
same reference. 

The revision of the scale resulted in marked 
improvement. The factor relating the two scales 
appeared to be 10-'’, but certain difficulties 
remain. The recent paper by Smallwood,‘ in 
which a like correlation has been made, induces 
the author to point out his correlation and cer- 
tain ways in which it may be applied. 

Table I lists the values of the coordinates for 
the various elements on the two electronegativity 
scales. The A! values are taken from Pauling’s* 
paper. He has calculated these values from the 
increment A, of observed bond energy in excess 
of calculated bond energy. In the latter case, the 
additivity of the energies of normal covalent 
bonds has been assumed. 

The uw values have been calculated from the 
observed electric moments of the hydrogen 
compounds. The moment of the molecule is 
directly that of the bond in the case of the 


‘ Zeits. f. physik, Chemie B19, 242 (1932). 
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TABLE I. Coordinates of elements on the two electronegativity 











scales. 
A} u X10 A} uX 108 

H 0.00 0.00 Br 0.75 0.78 
As -- 0.10 Cl 0.94 1.03 
y 0.10 0.36 N 0.95 1.04 

I 0.40 0.38 O 1.40 1.32 
Cc 0.55 = F 2.00 — 

S 0.43 0.63 








hydrogen halides, and since hydrogen is taken 
as the point of reference the value of the moment 
fixes the position of the halogen on the electro- 
negativity scale. In the case of oxygen and sulfur 
the values can be calculated only with a knowl- 
edge of the bond angles. The 90° angle taken, 
places the O—H bond at 1.32 on the scale in 
good agreement with the value given by Pauling. 
The position of sulfur at 0.63 is not in accord 
with the A? value of 0.43. Pauling has accepted 
the tandem structure for S2Cle in determining 
the position of sulfur. There is evidence against 
this structure. Objection has also been raised 
against the tandem structure for hydrogen 
peroxide.® The value given by Pauling for the 
O-—O bond would be substantially the same, 
however, if obtained from other sources. 

In the case of the fifth-group elements con- 
siderable data are obtainable regarding the bond 
angle.? Somewhat arbitrarily the angles 100°, 
98°, 96° and 94° have been assumed by the 
author as the bond angles of N—H, P—H, 
As—H and Sb—H, respectively. 100° is probably 
too low for ammonia, but it is nearer the true 
value than the 90° angle assumed by Smallwood. 
The decrease in the angles as one proceeds to 
higher weights is in line with the views of 
Pauling.® 

The position of carbon will not be given at this 
time as it presents difficulties. The substituted 
halides of methane follow no discernible pattern. 
Whether the difficulties arise from special proper- 


5 Matossi and Aderhold, Zeits. f. Physik 68, 683-95 
(1931). 

6 Maass, Can. J. Res. 781 (1932). 

7™L. Pauling, J. Am. Chem. Soc. 53, 1380 (1931); Wierl, 
Ann. d. Physik 8, 521 (1931); Dennison and Uhlenbeck, 
Phys. Rev. 41, 313 (1932); D. S. Villars, Chem. Rev. 11, 
369 (1932). 
8 L. Pauling, J. Am. Chem. Soc. 53, 1380 (1931). 
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ties of carbon, or from the circumstance that in 
every polar compound of carbon involving bonds 
of this type at least three species of atoms are 
included, is not known, but the latter appears 
the more probable. To use varying bond angles 
or induced polarization as an explanation is to 
defeat the purpose of this report. The electric 
moment values of the mixed halides of the fifth- 
group elements would furnish highly desirable 
information in this connection. 

The consistency of the electronegativity scale 
is indicated by the following moments of mole- 
cules calculated from bond moments taken from 
Table I and the bond angles previously used. The 
moment of the Sb—H bond has been assumed to 
be zero in calculating values for the antimony 
halides. 

From several viewpoints, it would be desirable 
to have electric moment data on the compounds 
resulting from the reaction of the halides with 
each other as these offer values for the group 
moments independent of any bond angle. The 
uX10'8 values calculated for the halides of the 
fifth group and recorded in Table III also possess 


TABLE II. Comparison of calculated and measured electric 











moments. 
uX108 
Molecule Calculated Measured Authority 
SCl, 0.51 0.56 J. W.S? 
PCl; 0.98 0.8 B. & E."° 
0.9 1.W.5* 
PBr; 0.62 0.61 B. & E. 
AsCl; 1.44 1.97 B. & E. 
AsBr; 1.05 1.66 B. & E. 
SbCl; 1.72 3.1 F. & M.! 
SbBr; 1.31 2.4 F. & M. 
SbI; 0.62 0.4 See refer 
ence 13 








9]. W. Smith, Proc. Roy. Soc. A138, 154 (1932). 

10 Bergmann and Engels, Zeits. f. physik. Chemie B13, 
232 (1931). 

11 J. W. Smith, Proc. Roy. Soc. A136, 256 (1932). 

2A. L. Ferguson and M. G. Malone, unpublished 
measurements. 

18 The value of 0.4. 10-8 is listed by Smyth (Dielectric 
Constant and Molecular Structure) for SbI3;. He attributes 
the value to Werner, Zeits. f. physik. Chemie B4, 371 
(1929). Werner did not measure this compound but quotes 
the value from Williams, Phys. Zeits. 29, 684 (1928). 
Williams and Allgeier, J. Am. Chem. Soc. 49, 2416 (1927) 
made an unsuccessful attempt to measure SbI; in benzene. 
The value 0.4 X 108 is doubtless too low. 
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TABLE III. »X10!8 values for the halides of the fifth-group 











elements. 
Sb As P N 
F; 4 = i ~~ 
Cl; 1.7 1.44 0.98 0.0 
(3.1) (2.0) (0.85) 
Br3 1.3 1.05 0.62 0.36 
(2.4) (1.66) (0.6) 
I; 0.62 0.43 0. 0.92 


(0.4) 








considerable interest. Where measured values 
are known they are given in parentheses just 
below the calculated values. . 
In view of the discrepancy between calculated 
and measured values, the predicted moments 
may be given only a semi-quantitative signif- 
icance. That the minimum moment of the iodides 
is indicated for the phosphorous compound and 
not the nitrogen, as hasty consideration might 


have indicated, illustrates a use of the electro- 
negativity scale. Similarly, the moment of the 
halogen compounds with nitrogen will increase 
in the order chloride, bromide, iodide instead of 
decreasing as occurs in the other columns. 

The perfection of two independent measures 
of electronegativity, one based on thermal data 
and one on dipole moment data, would be of 
considerable value. Not only would one supple- 
ment the other where experimental values were 
incomplete, but from the complete scales bond 
angles in polyatomic atoms might be determined. 
Before this can be done a clearer definition of 
electronegativity must be given. 

The author expresses his appreciation to A. L. 
Ferguson and M. G. Malone for the use of data 
not yet published in the usual way. Dr. Ferguson 
has also encouraged the writer to present this 
work at this time. 








MARCH, 1933 


JOURNAL OF CHEMICAL PHYSICS 


VOLUME 1 


The Dipole Moments and Structures of Certain Long-Chain Molecules 


C. P. SMytH AnD W. S. WALLS, Princeton University 
(Received December 30, 1932) 


The dipole moment of isoprene is measured and found to 
differ so little from zero that it cannot play a part in the 
tendency of the substance to polymerize. The moments of 
hexamethylene and nonamethylene bromide are deter- 
mined and that of trimethylene bromide is redetermined. 
Consideration of these values with those previously 
measured for other long-chain dibromides indicates that 
the moment of trimethylene bromide may be increased by 


the existence of a van der Waals attractive force between 
the two bromine atoms. The orienting effect of the two 
C—Br dipoles upon one another becomes negligibly small 
when the carbon chain between them is lengthened from 
six to nine atoms. The results are best explained by the 
supposition that the carbon chain is an extended zigzag 
structure which is not infrequently bent by rotation 
around one or more C—C bonds. 





INCE the completion of two series of meas- 

urements of the dipole moments of poly- 
methylene bromides,! Professor C. S. Marvel of 
the University of Illinois has very kindly lent us 
samples of hexamethylene bromide and non- 
amethylene bromide, which are here used to fill 
the two most important gaps left in the data on 
the polymethylene bromides. Because of the 
importance of isoprene in the building up of 
complex molecules, it seems not inappropriate 
to include measurements on it with those on the 
two long-chain molecules. Through the kindness 
of Dr. Thomas Midgley, Jr., and of Dr. Albert 
L. Henne, a sample of isoprene prepared by the 
latter was obtained for the measurements. As a 
check upon the earlier work, trimethylene 
bromide was measured in heptane solution. 


PURIFICATION OF MATERIALS 


Heptane 
Heptane was obtained as in earlier work.” 


Isoprene 

Isoprene was obtained from Dr. Albert L. 
Henne of Ohio State University as a mixture of 
60 percent isoprene with amylenes. For its 
molecular weight a weighted mean, 68.87, was 
used. mp* 1.40995. 


1C. P. Smyth and S. E. Kamerling, J. Am. Chem. Soc. 
53, 2988 (1931); C. P. Smyth and W. S. Walls, ibid. 54, 
2261 (1932). 

2C. P. Smyth and W. S. Walls, J. Am. Chem. Soc. 54, 
1854 (1932). 


Trimethylene bromide 


The material used was some of that purified ' 


for earlier work,' possessing the same refractive 
index as at the time of the earlier measurements. 


Hexamethylene bromide 

A somewhat impure sample, b.p. 92-93° (3 
mm) kindly lent by Professor C. S. Marvel of 
the University of Illinois was twice fractionally 
distilled under reduced pressure. After this 
material had been used for solutions the solvent 
was distilled off and the residual hexamethylene 
bromide was twice fractionally distilled under 
reduced pressure. A small sample of the middle 
fraction was used for making additional solu- 
tions, the polarizations of which checked those 
obtained for the first sample. 


Nonamethylene bromide 

A pure sample, b.p. 127-129° (7 mm) ob- 
tained from Professor Marvel was fractionally 
distilled slowly under reduced pressure, the 
middle fraction distilling at 103-105° (2 mm) 
being collected; mp”® 1.49521. 


EXPERIMENTAL RESULTS 


The densities and dielectric constants of the 
solutions of the polar substances in the nonpolar 
were measured with the apparatus and methods 
previously employed,? a wave-length of 1000 
meters being used in the dielectric constant de- 


°C, P. Smyth and W. S. Walls, J. Am. Chem. Soc. 53, 
527 (1931). 
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terminations. In Table I the first column gives 
the mole fraction c, of the polar substance in the 
solution and the succeeding columns give the val- 
ues at 25° and 50° of the dielectric constants e, the 
densities d of the solutions and the polarizations 
P. of the polar substance. The values given 
opposite c.=0 are those for the solvent. P2 has 
been extrapolated to c.=0 in the usual manner 
to obtain P.,, from which the moment is calcu- 
lated as u= 0.0127 X 10-8 (P,,— MRp)T }'. MRop, 
the molar refraction for the D sodium line, was 
determined directly for isoprene and trimethylene 
bromide and, for the two long-chain bromides, 
was calculated as the sum of the atomic refrac- 
tions. The values of these quantities are given 
in Table II, and, in Table III, the moments 
obtained for these dibromides are listed with 
those previously obtained by Smyth and 


TABLE I. Dielectric constants, densities and polarizations. 








Heptane-isoprene (60%) +amylenes (40%) 
(Mean M=68.87). 








€ d P, 

C2 = = rs 
0.00000 1.920 0.6795 (34.58 = P,) 
0.32132 1.959 .6770 27.27 
0.44984 1.978 .6760 27.29 
0.73520 2.028 .6737 27.30 
1.00000 2.087 .6715 27.28 

Heptane-trimethylene bromide. 
€ d P, 
Co 25° 50° ro 50° 2° 50° 





0.00000 1.920 1.883 0.6795 0.6577 (34.58 34.62=P,) 
03248 2.025 1.978 .7089 .6865 114.0 109.1 
07650 2.178 2.114 .7499 .7267 113.0 107.7 
13707 2.406 2.314 .8080 .7836 111.0 105.5 
18568 2.603 2.490 .8563 .8312 109.0 104.1 





Heptane-hexamethylene bromide. 





0.03148 2.065 2.014 0.7106 0.6884 162.0 155.0 
06780 2.239 2.174 .7467 .7239 158.9 153.6 
10946 2.447 2.365 .7877 .7643 155.9 151.5 
-13438 2.579 2.485 .8123 .7886 154.8 150.2 





(Middle fraction of recovered hexamethylene bromide.) 





0.02898 2.054 2.004 0.7082 0.6860 162.6 155.4 
05828 2.195 2.130 .7371 .7144 160.9 153.7 
05841 2.196 2.131 .7374 .7147 160.8 153.8 





Heptane-nonamethylene bromide. 





0.02368 2.041 1.994 0.7033 0.6816 191.3 184.5 
04631 2.159 2.101 .7255 .7035 189.8 183.2 
06669 2.262 2.195 .7452 .7229 186.5 180.4 
09153 2.393 2.313 .7689 .7464 184.6 178.4 
09891 2.429 2.346 .7757 .7532 183.3 177.3 














TABLE II. Refractions, orientation polarizations and electric 
moments. 








P,.—MRp pu X10" 


= 2 TF 





Isoprene (60%) +amy- 


lenes (40%) 25.41 1.87 (0.3) 
Trimethylene bromide 31.2 84.8 79.0 2.02 2.03 
Hexamethylene bromide 45.3 119.0 111.4 2.39 2.41 
Nonamethylene bromide 59.3. 135.2 128.2 2.55 2.58 








TABLE III. Dipole moments of polymethylene bromides. 











uX10"8 

Compound Formula Solvent 25° 50° 
Ethylene bromide Br(CH2)2Br Heptane-1.02 1.05 
Ethylene bromide Benzene 1.50 1.52 
Trimethylene bromide Br(CH2);Br Heptane 2.02 2.03 
Trimethylene bromide Benzene 1.97 1.98 
Tetramethylene bromide Br(CH2)sBr Heptane 1.96 2.01 
Tetramethylene bromide Benzene 2.00 2.03 
Pentamethylene bromide Br(CH:2);Br Benzene 2.25 2.28 
Hexamethylene bromide Br(CH:2)«Br Heptane 2.39 2.41 
Nonamethylene bromide Br(CH2),Br Heptane 2.55 2.58 
Decamethylene bromide Br(CH2)iwBr Benzene 2.54 2.56 








Kamerling for ethylene bromide and by Smyth 
and Walls for the higher dibromides. 


DISCUSSION OF RESULTS 


Although the formula of isoprene, CH,=CH 
-C(CH3) = CHe, gave reason to expect a moment 
too small to detect by means of measurements 
on the liquid, it seemed desirable to establish 
beyond all doubt that its ease of polymerization 
was not caused by the moment of the molecule 
as a whole. The apparent mean moment, 
0.310-'8, found for the mixture of isoprene 
and amylenes is so small that it may arise 
entirely from neglect of the atomic polarization. 
The fact that propylene* and a-butylene® have 
moments of 0.37 x 10~!8, while ethylene® has zero 
moment, shows that the double bond itself has 
no dipole moment, although it may give rise to 
a moment if unsymmetrically placed in a 
molecule. The maximum moment to be expected 
of an amylene is evidently about 0.4 10-'§ and 
the mixture would, doubtless, have a decidedly 
smaller moment. The experimental result shows 


*K. B. McAlpine and C. P. Smyth, unpublished 
measurements, 

5C. P. Smyth and C. T. Zahn, J. Am. Chem. Soc. 47, 
2501 (1925). 





202 Cu Ps 


a mean moment for the isoprene-amylene mix- 
ture between 0 and 0.310-!8 and, therefore, a 
moment for isoprene lying between approxi- 
mately the same limits, as would be predicted 
from the fact that the small moments arising 
from the double bonds at each end of the mole- 
cule would partially cancel one another. It is 
evident, therefore, that the dipole moment of the 
molecule is too small to play any appreciable 
part in the tendency toward polymerization. 

The moment of a polymethylene bromide 
molecule may, as a fair approximation, be 
regarded as the resultant of the moments in the 
two C—Br lines and is thus dependent upon the 
mutual orientation of these bonds. The difference 
between the moment of ethylene bromide in 
heptane and in benzene has been attributed to 
the weakening by the benzene of the forces 
between the two C—Br dipoles, which, by lessen- 
ing the potential energy between the dipoles, 
permits of their turning about the C—C axis 
farther out of the trans position.® 

The large energy due to the close proximity of 
the bromine atoms in the small molecule tends 
to keep the dipoles near this trans position, thus 
causing the moment of ethylene bromide to be 
much lower than the moments of the longer 
dibromides. Since, in the cases of trimethylene 
and tetramethylene bromide, the moments are 
almost the same in the two solvents, it appears 
safe to neglect the difference of solvent in the 
case of the longer molecules, where the effect of 
solvent should be even less. It was shown in 
the earlier work that, with one exception, the 
moments could be satisfactorily explained by 
the assumption of freedom of rotation about all 
the C—C bonds to give a great variety of forms 
to the molecules, those forms in which the dipoles 
point more or less in the same direction being 
less probable in the shorter molecules and the 
moments, consequently, being smaller than the 
calculated values, the difference between ob- 
served and calculated values lessening as in- 
creasing length of the chain lessens the mutual 
potential energy of the two dipoles. However, 
the closeness of the moment of tetramethylene 
bromide to that of trimethylene bromide was 


6C, P. Smyth, R. W. Dornte and E. Bright Wilson, Jr., 
J. Am. Chem. Soc. 53, 4242 (1931). 





SMYTH AND W. S. 





WALLS 


not accounted for by this explanation. It was 
pointed out that, if the carbon chain retained a 
rigid zigzag structure with freedom of rotation 
only about the end C—C bonds, the molecules 
with an even number of carbon atoms in the 
chain should have a moment 0.1X107'® lower 
than those with an odd number of carbons. This 
would reduce the rise in moment from the three- 
carbon to the four-carbon compound caused by 
increasing separation of the two dipoles and 
thus keep the values close together. In similar 
fashion, the moments of the five- and six-carbon 
compounds should be close together, but the 
moment of hexamethylene bromide, now avail- 
able to test the hypothesis, is decidedly higher 
than that of pentamethylene bromide. Because 
of its greater length, the six-carbon chain offers 
greater opportunity for some departure from 
the rigid zigzag structure than does the four- 
carbon chain, but the marked difference between 
the moments of the five- and six-carbon com- 
pounds renders improbable an explanation of the 
closeness of the three- and four-carbon values in 
terms of an alternating effect caused by a rigid 
zigzag chain. 

In our earlier paper, the mutual potential 
energies of the two dipoles in the different 
molecules were considered, but discussion of the 
energies due to other interatomic forces was 
postponed. Eyring’ has treated such forces in 
calculating the dependence of the potential 
energy of the ethane molecule on rotation about 
the C—C bond. In the present case, in addition 
to the forces between the dipoles, which affect 
their positions relative to one another in such a 
way as to reduce the resultant moment, there is 
a repulsive potential, such as Re~6”, which 
vanishes for distances R appreciably greater 
than kinetic theory diameters and a van der 
Waals potential, —C/R®, which gives rise to 
attraction. 

The constant C for an atom is easily calculated 
by the method of London.’ In trimethylene 
bromide, rotation around the two C—C lines, 
which make an angle of 110° with each other, 
would cause the distance between the bromine 
nuclei to vary from about 6A down to that of the 


7H. Eyring, J. Am. Chem. Soc. 54, 3191 (1932). 
’ F, London, Zeits. f. physik. Chemie (B) 11, 222 (1930). 
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closest possible approach of the two atoms to 
each other. If the two atoms lay in the same 
plane on the inside of the C—C—C angle, the 
nuclear separation of the two bromines would 
be only 1.53A. As this distance is considerably 
smaller than the diameter of a bromine atom,’ 
it is evident that this position and positions 
close to it are impossible in the trimethylene 
bromide molecule. As these are positions in 
which the two dipoles are not far from pointing in 
the same direction and thus giving a large 
moment to the molecule as a whole, it is evident 
that the moment observed for the molecule will 
be considerably smaller than that calculated on 
the assumption that all positions of rotation 
around the two C—C lines are equally probable. 
On the other hand, when the internuclear sepa- 
ration is a little greater than the kinetic theory 
diameter for the bromine atom, the van der 
Waals potential will become more important than 
the repulsive potential and there will be an 
attractive force between the bromine atoms, 
tending to overcome the repulsion between the 
dipoles pointing more or less in the same direc- 
tion. There should, therefore, be a value of the 
internuclear separation at which the attractive 
and repulsive forces should approximately com- 
pensate one another. The dipoles would then 
oscillate about this position of minimum poten- 
tial energy rather than about that in which the 
two point in opposite directions to give zero 
moment to the molecule and a minimum dipole 
potential energy. 

The compensation between the repulsive and 
attractive forces of the bromine atom in ‘the 
absence of dipole forces is shown by liquid 
bromine, in which calculation from the density 
shows that the average volume occupied by a 
molecule at 0° is equal to that of a cube with an 
edge of 4.35A. In terms of the dimensions given 
by Pauling, the Bre molecule would be roughly 
a cylinder of diameter 2.28A with hemispherical 
ends and total length 4.56A. These dimensions 
are, of course, decidedly smaller than those given 
by kinetic theory, from which Smyth and 
Rogers!” calculated a diameter 3.2A for the 
bromine atom. This would extend the hypo- 





*L. Pauling, Proc. Nat. Acad. Sci. 18, 293 (1932). 
“C. P. Smyth and H. E. Rogers, J. Am. Chem. Soc. 52, 
2227 (1930). 


thetical surface of the molecule calculated from 
the Pauling dimensions to a distance of 1.6A 
from the line joining the two bromine nuclei. 
If the liquid is regarded as made up of molecules 
of this shape packed parallel to one another in 
layers, it may be calculated that the average 
separation of the nearest bromine nuclei in 
adjacent molecules is about 3.6A. This very 
rough figure indicates a distance just above 
which the van der Waals attractive forces should 
tend to compensate the dipole repulsive forces 
in trimethylene bromide and give a point of 
minimum energy for the dipoles. With the bro- 
mines at such distances apart in the molecule, 
the dipoles make such angles with one another 
as to give considerable resultant moments to the 
molecules. This should cause the observed 
moment to be nearer that calculated on the 
basis of equal probability for all positions of 
rotation than would be expected in view of the 
shortness of the carbon chain. In other words, it 
is not that the tetramethylene bromide moment 
is low as previously supposed, but that the tri- 
methylene bromide moment is unexpectedly 
high. 

Although the van der Waals attraction may 
play a part in the tetramethylene bromide 
molecule, it must be a small part unless the 
carbon chain is distorted enough by rotation 
around the C—C bonds to bring the bromines 
close together. It seems probable that only in 
trimethylene bromide are the distances between 
the bromines commonly of such a size as to 
make the van der Waals attraction between 
them important. It is also of interest to note 
that, if the two dipoles in trimethylene bromide 
are oscillating about a position in which the 
resultant of their two moments is considerable, 
increased oscillation accompanying rising tem- 
perature will not greatly affect the observed 
moment. Actually, the values in Table III show 
a negligible increase in moment from 25° to 50°, 
and the earlier and less accurate determinations 
by Smyth and Kamerling showed an apparent 
increase in moment with temperature decidedly 
smaller for trimethylene bromide than for penta- 
methylene. 

Since the van der Waals potential varies 
inversely as the sixth power of the distance, it is 
negligible in the molecules containing more than 
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four carbons unless the ends of the chain are 
brought close together through rotation around 
the C—C lines. However, the potential energy 
between the two dipoles, which varies inversely 
as the cube of the distance, is still not incon- 
siderable when they are separated by a rigid 
zigzag chain of five carbon atoms. In connection 
with our earlier paper, it was calculated that the 
potential energy between the two dipoles when 
both were pointing upward in a sort of cis posi- 
tion was 1X10-“ erg greater than when one 
pointed downward in the opposite direction as 
in a trans position. As this difference is approxi- 
mately 1/4k7, it was evident that the trans 
position of zero moment was the more probable 
and that the observed moment would still be 
lower than the calculated and would show some 
increase with temperature as observed by Smyth 
and Kamerling. As the separation of the dipoles 
was approximately doubled in decamethylene 
bromide, it was evident that the difference 
between the potential energies of the cis and 
trans positions was about 1/8 of that in penta- 
methylene bromide, or about 1/32 kT, an energy 
too small to have any appreciable orienting effect 
upon the dipoles. In conformity with this, the 
observed moment was found to differ by less 
than the uncertainty of the calculation from the 
value calculated on the assumption of the equal 
probability of all positions of rotation, and was 
further found to show little or no variation with 
temperature. 

The approximate validity of the potential 
energy calculation receives striking confirmation 
from the additional moments in Table III. The 
moment of hexamethylene bromide is higher 
than that of pentamethylene because the dipoles 
are farther apart and have, therefore, a lower 
mutual potential energy. However, the difference 
in potential energy between the cis and trans 
positions is still not wholly negligible in com- 
parison to kT and is still, therefore, sufficient 
to make a little less probable the cis positions, 
which have greater moment, the resultant 
moment of the molecule, consequently, being 


SMYTH AND W. S. 


WALLS 





somewhat below the calculated value. In non- 
amethylene bromide, however, the difference in 
potential energy between the cis and trans 
positions is so small in comparison to kT as to 
produce no lowering of the moment, which is 
practically identical with that of decamethylene 
bromide and in satisfactory agreement with the 
calculated value. In the calculation of the result- 
ant moment of the two dipoles, a moment of 
1.9X10-'8 was assigned to each dipole and 
located in the C—Br line 1.7A from the carbon 
nucleus. The location of the moments is a rough 
approximation which makes approximate the 
calculation of the potential energy between 
them. In hexamethylene bromide the difference 
in potential energy between the cis and the trans 
position is about 1/7 kT, which is evidently still 
sufficient to have a small orienting effect on the 
dipoles. The difference in nonamethylene bro- 
mide drops to about 1/23 kT, which, as would be 
expected, is evidently insufficient to produce 
detectable orientation. 

It must be concluded that the moments in 
Table III give no positive evidence of an alter- 
nating effect such as might be expected if the 
carbon chain were a rigid zigzag structure, 
although, in the earlier work, the seemingly low 
moment of tetramethylene bromide was regarded 
as an indication of such a structure. However, 
the small difference, 0.110-'8, calculated be- 
tween the molecules with an even-numbered 
carbon chain and those with an odd would be 
lessened by any bending of the chain and more 
or less obscured by the lowering of the moment 
because of the mutual potential energies of the 
dipoles and because of inductive effects. Al- 
though, as pointed out in the earlier work, the 
moments do not decide definitely between an 
extended zigzag carbon chain and one in which 
free rotation around the C—C bonds gives a 
great variety of forms to the chain, the results 
are best explained by the supposition of an 
extended zigzag chain which is not infrequently 
bent by rotation around one or more links in the 
chain. 
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The Limiting Laws of the Interionic Attraction Theory of Strong Electrolytes 


PIERRE VAN RYSSELBERGHE, Chemistry Department, Stanford University 
(Received December 3, 1933) 


The various methods of obtaining the limiting laws of 
strong electrolytes are briefly discussed. It is shown that the 
dependence of these limiting laws on the temperature, the 
dielectric constant of the solvent and the volume of the 
solution can be deduced from the Gibbs-Helmholtz equa- 


tion and the virial theorem of Clausius. The Debye- 
Hiickel and the Kramers limiting laws are found to be of 
the same type as these general limiting laws, with numerical 
coefficients obtained thanks to the use of the Poisson 
equation. 





INTRODUCTION 


HE common features of the various forms 

of the interionic attraction theory of ionized 
solutes are: (1) The assumption of complete ion- 
ization at low concentrations (sometimes cor- 
rected by the consideration of a partial associ- 
ation of ions of opposite signs)'; (2) the assump- 
tion that the departure from the properties of 
perfect solutions is entirely due to the Coulomb 
attractions and repulsions between ions. 

The first step in any theory of ionized solutes 
is either the direct calculation of the correction 
that has to be introduced in the expression for 
the free energy of the system on account of those 
electrostatic forces; or the calculation of the 
energy involved in the so-called “charging 
process” of Debye.? Either one of these calcu- 
lations presents difficulties of theoretical and 
mathematical nature. Milner® tried to solve the 
problem directly by calculating the virial of a 
mixture of ions but his calculations are extremely 
intricate and contain a certain number of approx- 
imations. The method of Debye and Hiickel is 
based upon the Poisson equation of electrostatics 
which is only strictly valid for a continuous 
distribution of electricity. This procedure is 
rendered uncertain and weak on account of 


*N. Bjerrum, Det kgl. danske Vidensk. selskab, Math.- 
fys. Meddelelser VII, No. 9 (1926); W. Nernst, Zeits. f. 
Elektrochemie 33, 428 (1927); Zeits. f. physik. Chemie 
135, 237 (1928) ete. 

*P. Debye and E. Hiickel, Phys. Zeits. 24, 185 (1923); 
25, 97 (1924); E. Hiickel, Ergebn. d. exakt. Naturwiss. 3, 
199 (1924); W. Orthmann, ibid. 6, 155 (1927). 

Y *S.R. Milner, Phil. Mag. 23, 551 (1912); 25, 743 (1913). 


fluctuations, a point discussed by Fowler who 
considers the Debye-Hiickel theory as ‘‘empiric- 
ally successful.’”” The more exact form of the 
Debye-Hiickel theory developed by Gronwall,® 
by Gronwall, La Mer and Sandved* and by 
La Mer, Gronwall and Grieff,’ although much 
more satisfactory than the original form of the 
Debye-Hiickel theory, also suffers from fluctu- 
ations of an unknown order of magnitude since 
it is also based on the Poisson equation, for 
which an exact solution is now obtained. The 
only parameter of the theory, the mean ionic 
diameter, is found to have more consistent values 
than in the original Debye-Hiickel theory, 
where absurdly small or even negative values 
were often found. 

A remarkable solution of the problem, ‘‘avoid- 
ing the pitfalls of Debye and Hiickel’s method 
and the tedium of Milner’s’’ (Fowler) has been 
obtained by Kramers.’ It consists of a direct 
calculation of the correction in the free energy 
due to the electrostatic forces, by means of 
Gibb’s phase integral. The Poisson equation is 
not used (at least not as the fundamental point 
of departure of the theory) and hence fluctuations 


*R. H. Fowler, Statistical Mechanics, pp. 194-197, 316- 
325, Cambridge, 1929. 

5 T. H. Gronwall, Proc. Nat. Acad. Sci. 13, 198 (1927); 
Annals of Math. 28, 355 (1927). 

* T. H. Gronwall, V. K. La Mer and K. Sandved, Phys. 
Zeits. 29, 358 (1928). 

7V. K. La Mer, T. H. Gronwall and L. J. Grieff, Phys. 
Chem. 35, 2245 (1931). 

8H. A. Kramers, Koninkl. Akad. van Wetenschappen te 
Amsterdam, Afdeeling Natuurkunde, XXXV, No. 10, 1153 
(1926). 
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are avoided. It is shown that up to definite 
concentrations depending on the valence type 
of the electrolyte the part of the free energy due 
to the electrostatic forces, the activity coefficient 
and the osmotic coefficient are independent of 
the size of the ions and the Debye-Hiickel 
limiting laws (under their simplest form, corre- 
sponding to negligibly small ionic radii) are 
obtained for the case of extremely low concen- 
trations. Kramers further shows that the treat- 
ment given by Debye and Hiickel for higher 
concentrations (finite size of the ions taken into 
account) cannot be significant because at those 
concentrations association takes place and rea- 
sonings based on the existence of atmospheres of 
free ions are invalid. It is to be noticed that the 
Kramers theory is developed without the use of 
any parameter. Towards the end of his reasoning, 
however, Kramers uses the Poisson equation in 
order to determine the numerical value of a 
constant. 

The object of this note is to show that it is 
possible to obtain directly the general form of the 
limiting laws of strong electrolytes from the 
combined use of the Gibbs-Helmholtz formula 
and the virial theorem of Clausius. It is then 
shown that the Debye-Hiickel and the Kramers 
limiting laws are particular forms of these general 
laws. 


LimITING LAWS OF PERFECT STRONG 
ELECTROLYTES 


We call “perfect strong electrolyte’ a com- 
pletely ionized solute satisfying the following 
conditions: (1) The ions are point charges. (2) 
The departure of the osmotic properties from 
those of an ideal solution is entirely and only due 
to the Coulomb forces between the ions. (3) The 
dielectric constant of the solvent is a function of 
temperature only. 

Let us call €;, €2, «++ €:, «++ the charges of the 
various ions, and yj, Wo, «++ Wi, °++ the electro- 
static potentials at the points occupied by the 
charges €1, €2, «++ €, **+ and due to all the other 
charges. Thus we have: 





1fe €2 €i-1 €i41 

vena [edt greek, OD 
Dir rai Ti-1, 4 i+1, ¢ 

Vii, Toiy *** Tei *** being the distances between 
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the various ions and the ion and D the dielectric 
constant of the medium. It is obvious that we 
are dealing with the time averages of the y,’s 
and ei S. 
The electrostatic energy U, of the system is 
given by 
Ue=2 Dili. (2) 


Let us call P the osmotic pressure of the solution, 
P the osmotic pressure of the corresponding 
ideal solution, A the Helmholtz free energy or 
maximum work of the system, A the same 
quantity for the ideal solution, V the volume of 
the solution. We have: 


P=—(0A/aV)r, (3) 
P=-—(dA/dV)r. (3’) 


A could be decomposed into A and W, the 
energy involved in Debye’s charging process’: 


A=A+4W. (4) 
From (3), (3’) and (4) we deduce: 
P=P-—(aW/aV)r. (5) 


The virial theorem of Clausius leads to the equa- 
tion of state: 


- 


PV=NRT+3>i> erecc'f (rei), (6) 


where N is the total number of particles in the 
volume V, k is Boltzmann’s constant, T is the 
absolute temperature, 7;; is the distance between 
two particles, 7 and k, f(ri:) is the force acting 
on the particle 7 on account of particle k. We 


have: _ 
P=NkT/V. (7) 


We deduce from (5), (6) and (7): 
—(0W/dV)r=(1/3V)-SsdXierecf(rei). (8) 
From (1) sad (2) we Gad: 
LeiLeereif (res) = (1/D) Didi eres(ever/rin’) 
~ LL wresf (res) = 2 Vise. (9) 
Eq. (8) gives then: 
—(AW/dV)r=(1/6V)- Die. — (10) 


But since W is the energy spent in charging the 
ions against the potentials y;, we have: 








In 


fol 
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w=: f ~ (11) 





Eq. (10) becomes then: 


0 «i 1 
|= f vdes |=—- Seed (12) 
0 6V 


av 
We have for the ion 7: 
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fi(ex/ V'/®) being an arbitrary function of ¢;/V'/®, 
According to Bjerrum’ the heat of dilution of 
an electrolyte is given by: 


AU=U.[1+(T/D)-dD/dT]. (16) 


AU is the heat evolved when the solution of 
volume V is diluted to an infinite volume. 

If the solution were perfect, AU would be zero. 
The Gibbs-Helmholtz formula gives: 














€i OY: 1 
-{ ial ee (13) AU=A-—T:(0A/aT)y. (17) 
0 OV 6V But 7 
Differentiating with regard to «;: H edited ieads (18) 
ence: 
vim —eeOYi/Oee—OV-OW/OV. = (14) | no A_T(@A/aT\y+W-T-(OW/aT)v. (19) 
) Integrating: —— 
> Wi= (1/e;) -file./ VV"), (15) A-—T(dA/dT)y=AU=0 (20) 
: we obtain 
T dD Ow 
) AU= u.[1+—-— |=w-7 (—) (21) 
D dT dT sy 
or: 
1S sea) [1+ —] . [va - rs) [x if . 2 
22 i€iWi oe a ae Pa a i Wide; |. ) 
) ; D dT 0 aT 0 
For the ion 7: 
T dD ci “OW: 
. Lew; | 1+—-— |= f bader-T dé. (23) 
D aT 0 0 OF 
Differentiating with regard to e¢;: 
q 
Ce (14- —) (19 43 ~~) 4 OY; 
> © <oneame Wi a = ae 24 
pare a wid 
4 or: 
0 T dD T dD} oy: OY: 
¥| 1 =a [14]. +2T-—. (25) 
D dT D dT J 0d¢; oT 
Integrating (D is a function of T only): 
yy i= T'D-}. f(e,/T'D!) (26) 
fo(e;/T*D#) being an arbitrary function of e;/T!D*. Let us expand expressions (15) and (26) into 
infinite series, according to Maclaurin’s formula. We get: 
1 &j 1 €? 1 e? 
v= —| HOF HOA HOF HO fees (27) 
) €; Vi 2! Vi 3! V3 
€; 1 «7 1 e? 
¥i:=TiD- | 1(0)+ fo’ (0) +—-—- f."(0) +—- fol""(0) ++: | (28) 
) T?D3 2! TD 3! TDi 
e One easily finds that f,(0) =0, fo(0) =0, f:/(0) =0. Hence, introducing the Boltzmann constant, so 











*N. Bjerrum, Zeits. f. physik. Chemie 119, 145 (1926); see W. Orthmann, reference 2, p. 181, 
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as to make the C’s dimensionless: 
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€i €? é; 
i=C +C2 Or prapcenmcpoemy medial 
DV (kRT)§DIV3 (kT)D?Vi 
or: (29) 
m=O €;” 
Wi= DUC 





m=” (RT) (m-DI2. Pmt 2. Yom41)/6 


The C,,’s are constants depending on the total 
number of ions, their valences, etc. Their deter- 
mination requires calculations of the type of 
those of Debye and Hiickel, of Milner and of 
Kramers. 

In the case of negligibly small ionic radii 
Debye and Hiickel found: 


Vi= —(e:/D)-((4re/DRT)- > Niz?2/V)}, (30) 


e being the elementary quantum of electricity, 
N; the total number of ions of the kind 7, 2; the 
valence of these ions. 

One sees that Eq. (29) reduces to (30) when: 


1=Cs=Cy=--- =(0 
and 
C2= —(4r-> Niz2/27)3. (31) 


The value of W obtained by Kramers is the 


following: 
W=3NkT In w(x). (32) 


N is the total number of ions, w(x) is a function 
whose dependence on x is defined as follows: 


w=e'(14+32)-, (33) 


(4r/9)x =P(1+32) (34) 
and 
x=[(2/DkRTN)-> Ni? }?- N/V. (35) 


When x is much smaller than 1, we have: 


w= 1—(2/3)(4rx)!—(5/18)-4rx—---, (36) 





—In w=(2/3)(4ex)!+34ex+-++, (37) 
W=-X(Ni2/3)-ex/D—--» (38) 


with x=((4re?/DkT)->-Nj22/V)*. The first 
term of W is the value found by Debye and 
Hiickel for this quantity in the limiting case of 
negligibly small ionic radii. 

Let us calculate y; from (32) and (37). We 
have: 


i= (2:/> Niz?)-OW/de. (39) 
This gives: 





y= —-—- 
D \DkT V 


€; 67e! (>> N27)? 


D (DkT)}? NV 


€; 4ré ane) 





(40) 


Comparing with Eq. (29) we see that: 


1=C3=Cy=Ce=--- =0 
and 
C2 = (47- > N,22/2,7)?, 
C;= —6r7- (> Niz,7)?/ Nz;4, etc. (41) 


From Eq. (29) we could deduce, if the C.,,’s 
were known, the osmotic coefficient g and the 
activity coefficient f by means of the following 
relations: 


g=P/P=1-(dW/aV)r, (42) 


(43) 





w=¥,|N e 


‘om (m+1)(RT)™-D/2- DimtD/2, Yomi) /6 


In f;=(1/kT)-OW/ONi. (44) 


The activity coefficient of the salt 7 will be 
given by: 


f= Lifer fri}”, (45) 


v; and v; being the numbers of ions 7 and 7 into 
which the salt 77 dissociates, v being the sum of 








Ci, meh 
, 


v; and v;. It is seen that g and f are given by 
series. If these series are reduced to their first 
terms, g and f will be given by expressions 
analogous to those of Debye and Hiickel, exactly 
as in the Kramers theory. This approximation is 
justified at very low concentrations and for not 
too small values of D. It happens that these first 
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terms are also obtained when the Poisson equa- 
tion is integrated according to the method of 
Debye and Hiickel (hyperbolic sine replaced by 
its argument), i.e., when in the case of 1—1 
electrolytes, for instance, the equation 


Ay = — (4ane/D) (e-*¥/*7 — vk?) (46) 
is arbitrarily replaced by 
Ay = 8rney/DkT. (47) 


This result cannot be considered as anything 
more than a coincidence, since La Mer, Gronwall 
and Sandved* have shown that if the exact 
integration of the Eq. (46) is carried out, the 
value of ¥; and hence that of the electrostatic 
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energy of the mixture of ions become infinite 
when the radius of the ions becomes infinitely 
small. This would tend to show that the Poisson 
equation as applied to electrolytes is meaningless 
in this case. 

Summing up we may say that the reasonings 
developed in this paper show how far it is 
possible to go in the solution of the problem of 
electrolytes without using more statistics than 
what is involved in the virial theorem and with- 
out using the Poisson equation. These reasonings 
also show that, insofar as their dependence on 
the temperature, the dielectric constant of the 
solvent and the volume of the solution is con- 
cerned, the limiting laws of Debye-Hiickel and 
of Kramers are essentially correct. 
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Wave Functions for the Ground State of Lithium and Three-Electron Ions 


E. BriGut WILSON, JR., Gates Chemical Laboratory, California Institute of Technology 
(Received January 16, 1933) 


Improved wave functions and ionization potentials have 
been obtained for the configuration 1s? 2s, by using the 
variation method with a variation function containing four 
parameters. The wave function for the whole atom is of the 
determinant form, built up of K and L single-electron 
functions. The K function used is hydrogen-like with one 
parameter, the effective nuclear charge, and the LZ function 
is of the form are~"*—e*, A simple rule is obtained by 


means of which it is possible to write down the wave 
function for any three-electron ion without further calcu- 
lation. The deviation of the calculated ionization potentials 
from the observed values is 0.9 percent for Li I and de- 
creases to 0.3 percent for C IV. Comparison is made with 
the results of other investigators and curves showing the 
wave functions and electron density functions are given. 





INTRODUCTION 


HE solution of the Schrédinger equation for 

the hydrogen atom and for the ground 
state of the helium atom! has been carried out to 
a high degree of exactness, but the best approxi- 
mation which has been obtained for the lowest 
state of lithium is not nearly so accurate. While 
it is probably not feasible at present to obtain a 
result for the energy which is superior to the 
experimental value in accuracy, it was thought 
worth while to test out on lithium various 
possible approximate solutions, partly because 
the better agreement with the experimental 
energy and the correspondingly improved wave 
function and electron-density function were con- 
sidered important in themselves and partly 
because this study would provide a basis for 
similar investigations of heavier atoms. 


GENERAL METHOD 


Slater? has shown how to build up a properly 
antisymmetric wave function for a many- 
electron atom by using a determinant whose 
elements are functions of the coordinates of a 
single electron only. In particular, if the inter- 
action of the electrons with each other is regarded 
as a perturbation and neglected in obtaining the 
zero“ order approximation to the solution of the 
wave equation, there results for lithium: 


1E. A. Hylleraas, Zeits. f. Physik 65, 209 (1930). 
2J. C. Slater, Phys. Rev. 34, 1293 (1929). 


210 


, A, Ao As 
vad, A: As (1) 
B, Bo Bs 


where A, is the 1s hydrogen-like wave function 
with positive spin for the first electron; A: is the 
same function with negative spin for the second 
electron, and B; is the 2s hydrogen-like function 
with plus spin for the third electron. 

The method adopted in this work was to seek a 
better approximation by the use of new functions 
for B based on the hydrogen-like functions but 
generalized by the introduction of parameters 
whose values were then determined by the 
application of the variation method. 

The solution of the variation problem: 


bE = 5 Sy*Hydr/ Sf y*ydr ]=0, (2) 


where y is the function varied and H is the 
Hamiltonian operator for the system, is equiva- 
lent to solving the wave equation 


Hy=Wy (3) 


with the usual boundary conditions. In particular 
if this variation problem is not completely 
solved but instead a trial function y used which 
contains a number of parameters which are 
varied until E is a minimum, then it can be 
shown’ that E is an upper limit to the energy, 
and y is assumed to be an approximate solution. 
The method which was used in this problem is 


*C. Eckart, Phys. Rev. 36, 878 (1930). 
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applicable to the lowest state of any multiplicity 
or total angular momentum. 


PREvious WorK ON LITHIUM 


Hargreaves‘ has applied the method of Hartree 
to lithium, but his results have not appeared in a 
very usable form. Hartree’s method of the 
“self-consistent field’”’ is not based on the vari- 
ation principle, but Slater® has shown the relation 
between these two procedures and in addition has 
pointed out that several important corrections 
are needed before the energy values obtained can 
be compared with the experimental ionization 
potentials. The curve Hargreaves publishes for 
the wave function cannot be used for comparison 
since he gives no scales and has not made the L 
function orthogonal to the K function, so that its 
shape is of little significance. However, he does 
give a table of the self-consistent field, and these 
values will be compared with other results in 
Table ITI. 

Eckart,’ using essentially the same method as 
is described in this paper, employed hydrogen- 
like functions with two parameters, the effective 
nuclear charges for the K- and L-shells. Thus if 
we use a, b, c, «++, to represent single electron- 
orbit functions without the spin factor and 
A, B,C, ---;A, B, C, ---, for the same functions 
with plus and minus spin, respectively, Eckart’s 
solutions are 


(K-shell functions) a =e-*’ 


4 
(L-shell functions) b=e—"’(nr—1). (4) 


There appears to be an error in his energy 
value for lithium so this has been recalculated 
and will be found in Table I. 

Guillemin and Zener* introduced an additional 
parameter (a) into b, obtaining thereby a con- 
siderable improvement in the energy. Their 
function is: 

b=e—(ar—1) (5) 


and values for the energy will also be found in 
Table I and of the parameters in Table II. In 
addition these authors varied the exponent of r 





‘J. Hargreaves, Proc. Camb. Phil. Soc. 25, 75 (1929). 
*J. C. Slater, Phys. Rev. 35, 210 (1930). 
' *V. Guillemin and C. Zener, Zeits. f. Physik 61, 199 
1930), 
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TABLE I. Energies and ionization potentials (atomic 
units: e*/do). 








Toniza- 
tion 
Total poten- 
energy Diff. % © tial Diff. % 





Lil 


Experimental —7.4837 — — 0.1983 — — 

4-Parameter —7.4192 0.0645 0.86 .1965 0.0018 0.91 

Guillemin- —7.4183 .0654 .87 .1956 .0027 1.36 
Zener 

Slater —7.4179 .0658 .88 .1953 .0030 1.51 

Hydrogen-like —7.3922 .0915 1.22 .1696 .0287 14.5 














Be II 
Experimental —14.3422 — — .6704 — — 
4-Parameter —14.2639 .0783 .55 .6663 .0041  .61 
Slater — 14.2584 .0838 .58 .6607 .0097 1.45 
B III 











Experimental — 23.476 — — 1.395 — — 

4-Parameter —23.363 .113 .48 1.390 .005 .36 

Slater —23.350 .126 .54 1.378 017 1.22 
CIV 

Experimental —34.778 —  — 2.3722 


4-Parameter —34.713 .065 .19 2.3650 .0072 .30 


Guillemin- — 34.698* .080 .23 2.3496* .0226 .95 
Zener 
Slater —34.690 .088 .25 2.3422 .0300 1.26 








*This is probably not quite the best value obtainable 
with this function. 


TABLE II. Parameter values, except for four-parameter 




















function. 
g n a é n a 

Lil Be II 
Guillemin-Zener® 2.688 0.630 5.56 3.688 1.158 5.88 
Slater 2.688 .630 — 3.682 109 — 
Hydrogen-like* 2.686 .888 — 3.70 142 — 

B III CIV 
Guillemin-Zener® 4.688 1.671 5.56 5.688 2.179 5.88 
Slater 4.676 1.52 — 5.672 1.95 — 
Hydrogen-like® 4.72 196 — 5.71 248 — 








in the coefficient, but found in this case that its 
best value was one, the hydrogen-like value. 
Zener’ and Slater® noted that a is fairly large 
and concluded that a good result would be 
obtained with the very simple function 


b=re" (6) 


7™C. Zener, Phys. Rev. 36, 51 (1930). 
8 J. C. Slater, Phys. Rev. 36, 57 (1930). 


This function was tried and, as is seen in Table I, 
found to be remarkably good considering its 
simplicity. It is better than the hydrogen-like 
function with screening constant used by Eckart, 
although, of course, not as good as Guillemin and 
Zener’s function, of which it is a special case. 


FUNCTIONS INVESTIGATED 


The type of function which has been used in 
this paper, a natural generalization of Guillemin 
and Zener’s function, is 


b=are" —et", (7) 


The same type has been used by Slater® in seeking 
analytical expressions for Hartree’s graphical 
wave functions. Results which are obtained by 
using (7) are listed in Table I under the heading 
“‘4-parameter’’ function. 

The following form was also tried but found 
not to give as good results as (7): 


b=are-"+re*". (8) 


Since the variation principle is applicable to 
the lowest P state of the atom, the following 
function was tried for the configuration 1s? 2p of 
lithium: 

c=(ye*"+e-")r cos 8. (9) 


However, the surprising result was obtained, on 
varying the parameters, that this function re- 
duced to the simple screening constant type used 
by Eckart. In other words, the best values of the 
parameters y, x and \ were such that y=0; «=X. 
Therefore, the calculated total energy for this 
configuration, —7.35039 e?/ao, is not different 
from Eckart’s value, but by combining this with 
the best calculated value for the ground state, 
—7.41915 e?/ao, a new result for the resonance 
potential of lithium is obtained, 0.06876 e?/ao, 
which is to be compared with the experimental 
value, 0.06794 e?/ao. The difference is 0.0082 e/a 
or 1.2 percent. 

The calculated energy of removal of the 2p 
electron from 1s? 2 of lithium is 2.2 percent 
greater than that for a 2p hydrogen atom, 
e?/8a9, whereas the experimental value is 4.2 
percent greater than e?/8a9. The discrepancy of 2 
percent is probably to be attributed to the effect 


9 J. C. Slater, Phys. Rev. 42, 33 (1932). 













































212 E. BRIGHT WILSON, JR. 


of polarization of the K-shell by the valence 
electron, the effect of which is not taken into 
consideration by any of the variation functions 
used in this paper, approximate calculations of 
the polarization energy” leading to somewhat 
larger values than the 2 percent needed. It is of 
interest that these calculations show, as sug- 
gested by Pauling and Goudsmit," that there is 
appreciable interaction of the valence electron 
with the core of the atom other than polarization 
even for ‘“‘nonpenetrating”’ orbits. 


DETAILS OF METHOD 


When the units a)=h?/(4r?ye”), 2Rhc=e?/ao 
and e are used for length, energy and charge, 
respectively, Schrédinger’s equation for the 
three-electron problem becomes 


Hy =(—3V°+V)y=Wy (10) 


where 
V=1/rio+1/re3+1/ri3—Z/r1—Z/re—Z/rs. (11) 


V? is the Laplacian for the coordinates of the 
three electrons, and Z is the atomic number of 
the atom or ion. As pointed out by Hylleraas,” 
application of Green’s theorem to the term in (2) 
involving the Laplacian yields for E: 


E=(M’'—2L’)/2N’ (12) 
where 
M'= Sf \grad y\*dr 
L'’=—-—fy*Viydr ¢. (13) 
N= Svtvdr 
Further, if the change in scale: 
v=tr; n=nlt; o=S/k (14) 


is made, it is possible to differentiate E with 
respect to £, solve for the minimizing value of & 
and insert this value in E. The result is 


E=-—L?/2MN (15) 


where L, M, N differ from L’, M’, N’ only in that 
the substitutions (14) have been made. 

The other parameters enter in such a compli- 
cated fashion that it is not feasible to obtain the 


101, Pauling and S. Goudsmit, The Structure of Line 
Spectra, McGraw-Hill, New York, 1930, p. 47. 
" Reference 10, p. 48. 
2 E. A. Hylleraas, Zeits. f. Physik 54, 347 (1929). 
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minimum by differentiation so resort was had to 
actually computing FE for various values of the 
parameters and graphically determining the 
minimum. This is a tedious process at best; the 
quickest and most certain method was found to 
involve constructing rough contour maps of the 
energy surface. An estimate was first made of the 
probable values of parameters. One of these 
(usually »’) was inserted in the formula for E 
and a survey of the ¢’, a’ plane made by giving ¢’ 
a number of fixed values and plotting E as a 
function of a’. Each of these curves then gives a 
section of the energy surface, with E thought of 
as a function of a’ and ¢’ only. From these a 
contour map can be constructed and an estimate 
made of the best values of a’ and ¢’ for the value 
of 7’ originally chosen. With these values of a’ 
and ¢’, 7’ is then varied until a minimum is 
found. A contour map may be constructed for the 
a’, ¢’ plane with this new value of 7’ and the 
best values redetermined. 


K-SHELL FUNCTIONS 


In all the cases discussed here, with the 
exception of Hargreave’s paper, the function 
used for the K electrons was the same as in 
Eckart’s work, i.e., a hydrogen-like function 
with one parameter, the effective nuclear charge, 
as in Eq. (4). It is well known that this is not the 
best function to represent the K electrons, some 
term such as one of those used by Hylleraas,” 
for example, to correct more exactly for the 
repulsion of the electrons giving a better result. 
However, the introduction of even the simplest of 
these interaction factors, such as (1+ /7r12), 
enormously complicates the problem in the case 


A, Az As 











| a - 1 
y= A, Az As=— 
63 63 

B, Bo Bs 


and in particular 8 may be so chosen that 
SA,(Bi+BA;)dr=0. (17) 


With this value of 8, the single-electron f unctions 
are mutually orthogonal, since A; and A, have 
Opposite spins. In addition, the functions a and 





A, A» A; 


Bi+BA, Be+BA2 B3+ BA; 
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of lithium and so was not regarded as practicable. 
Instead, the ionization potential was calculated 
by subtracting the energy of the Li* ion com- 
puted by using the same type of 1s wave 
function, the value of the energy so calculated 
being: 7.22266 e?/ay for Lit; 13.59766 e/a» for 
Bet+; 21.97266 e?/ao for B***; and 32.34766 e?/a 
for C++++, Although the total energy of the 
lithium atom is, therefore, in error by about 
0.065 e?/ao or 0.86 percent in all these examples 
due to the inaccuracy of the K functions, the 
assumption made that this error is nearly equal 
in the ion and in the atom is well justified by the 
close agreement between the experimental and 
theoretical ionization potentials as given in 
Table I. 


COMPARISON OF L FUNCTIONS 


The 2s functions which have been enumerated 
do not, at first glance, appear to be very similar, 
and it may seem strange that such different 
functions can all be even moderately good 
approximations to the true solution. The reason 
for the discrepancy is the use of single-electron 
wave functions for comparison instead of the 
wave function for the atom as a whole. Slater® 
has pointed out that it is possible to make the 2s 
function orthogonal to the 1s function by adding 
to the 2s a certain fraction of the 1s function. 
Since the wave function for the whole atom is a 
determinant, it possesses the property that the 
addition of the members of any row multiplied 
by any factor to the members of any other row 
does not change the value of the determinant. 
Thus: 


A, Ag As (16) 





b later tabulated have been normalized to 47, 
i.€., 
SAfdr= f BYdr=4r. (18) 


The orthogonalization and normalization make a 
and b unique (as long as functions of different 
spins are not combined), and it is thus possible to 
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compare the different approximations which have 
been obtained. 

As shown by Fig. 3, the four approximations 
already enumerated are remarkably near to each 
other. The only one showing any great deviation 
is the hydrogen-like function, which is also the 
one giving the poorest energy value. The simple, 
so-called ‘‘nodeless’’ function of Slater® is seen 
to be miscalled, since, when it is orthogonalized, 
it has a node in practically the same place as the 
more precise functions. 


RESULTS 


The parameter values obtained for the func- 
tions of Guillemin and Zener, Slater, and Eckart 
are found in Table II. In the case of the four- 
parameter function, the best parameter values 
fit into a very convenient rule if we introduce the 
idea of screening constants. The function for the 
L electron which has been used, Eq. (7), consists 
of two terms (before orthogonalization), and it is 
seen that the first term is the more important at 
large distances, because of the factor r and the 
smaller coefficient in the exponent, while the 
second term is the more important for small 
values of r. In the simple hydrogen-like function 
analogous to this, there is only one exponential, 
whose exponent is Zesr./2, or (Z —o)/2 if we define 
the quantity o, the screening constant. o thus 
represents the shielding effect of the inner 
electrons which reduces the attraction of the 
nucleus for the L electron. Likewise in our more 
complicated function the idea of the screening 
constant can be introduced, only here two 
different screening constants for the Z electron 
are needed, since there are two exponentials. The 
results which are obtained are that these screen- 
ing constants are essentially independent of the 
nuclear charge and, therefore, provide a con- 
venient method of summarizing the results of the 
variation problem. Furthermore, to a sufficiently 
close approximation, the screening constant for 
the inner part of the 2s function is zero, so that 
in a sense the four parameters originally used 
have been reduced to three. 

The parameter a can also be related for the 
ions of different atomic number. It is found that 
in all four cases investigated a’ comes so close 
to 0.5 that it is very convenient to assume that 








BRIGHT WILSON, JR. 



































~ 23302415 1 
| Kf 
-23.362500} “<i 
a _—_— 
~23.362525) 
~ 23.362550 
048 049 0.50 035i O52 


' 
Oo — 


Fic. 1. Curve showing the variation of E with the 
parameter a’ for B III, with the four-parameter function 
with o;=0.31, o2=1.67, ¢=Z/2. 


value for all three-electron ions. Fig. 1 shows the 
way E varies with a’ in the case of B III. 
Since a=a’t and §£=Z—«a,, if o; is the K-shell 
screening constant, therefore, 


a=0.5(Z —o,) =0.5(Z—0.31). (19) 


Introducing the values indicated above for the 
parameters 7, ¢, a and £, we obtain for the best 
wave function for the state 1s? 2s 2S of any three- 
electron atom or ion the expression 





v= (1/6!) p+ PA 1A oB;. (20) 
a= (§/r)'e& (21) 
b=k{ are” —e*"+ Bet} (22) 
where 
§=Z—-01 ¢=Z2Z/2 (23) 
n=(Z—o2)/2 a=0.5(Z—«;) 
= Bes _ Sax orthogonality 
ae ee 





oe 1 6 48a 
7 ff # (tty 
4808 168 4 
ae, 
(n+&)* (¢+6)? 


and Z= atomic number of ion; and o,=0.31, K- 
shell screening constant. ¢2= 1.67, L-shell screen- 
ing constant. The sum in Eq. (20) is over all per- 
mutations P of the numbers 1, 2, 3, odd permu- 
tations having the negative sign. (This is another 
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OF LITHIUM 
way of writing the determinant of Eq. (1).) The 
single-electron functions a and b as given above 
are individually normalized to unity and are 
mutually orthogonal. 

The energies obtained with the various func- 
tions are tabulated in Table I, which also con- 
tains the calculated ionization potentials. The 
experimental energy values are listed in the same 
table, together with deviations and percent 
deviations of the calculated and observed quanti- 
ties. It is seen that the poorest values are those 
obtained with the hydrogen-like function of Eq. 
(4), even the simpler function put forward by 
Slater and given in Eq. (6) being considerably 
better. The error in the ionization potential 
computed by the use of the four-parameter 
function is only from thirty to seventy percent of 
the corresponding error found when Guillemin 
and Zener’s function is used. 
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Fig. 2 is a plot of the four-parameter single- 
electron function for the 2s electron normalized 
to 4m and with 8 adjusted so that the Z and K 
functions are orthogonal. Fig. 3 gives 47r°b?, the 
distribution function integrated over all angles. 
On comparing the size of the two maxima in Li I 
and in hydrogen (2s state), the effect of the core 
electrons of lithium in reducing the probability of 
close approach to the nucleus by the valence 
electron is clearly seen, in agreement with the old 
concept of the penetrating orbit, in the inner 
segment of which the electron was speeded up by 
the increased effective nuclear charge, diminish- 
ing the time spent in traversing this segment. Fig. 
4 gives the total electron density function for 
Li I, integrated over all angles. The very steep 
rise at the K-shell and the more spread-out 
maximum of the single Z electron are shown. 
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Fic. 2. Plot of four-parameter orthogonalized single-electron wave function 6 for 2s electron, normalized to 4r. 


r in atomic units. Li I. 
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THE INTERNAL FIELD 


Hargreaves‘ has given a table of values for the 
self-consistent field for lithium, obtained graph- 
ically. This is expressed in terms of the “effective 
nuclear charge for field,’’ which is the number Z; 
such that the charge Ze placed at the position of 
the nucleus would give the same field (not 
potential) at the point in question as is actually 
found. In other words: 


Z;=Pr(dV/dr) (24) 


in atomic units. 
In Table III the values of Z; found by 
Hargreaves are listed. In addition there is also 
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given the field due to the nucleus and the two K 


TABLE III. Field due to K electrons. r in atomic units. 











r Z;—Hargreaves Z;—Eq. (24) Difference 
0 3 3 0 
0.1 2.963 2.965 —0.002 

on 2.807 2.811 — .004 

3 2.559 2.559 .000 

A 2.280 2.272 + .008 

6 1.773 1.748 + .025 

8 1.425 1.394 .031 
1.0 1.219 1.192 .027 
1.2 1.108 1.089 .019 
1.4 1.050 1.039 O11 
1.6 1.022 1.017 .005 
1.8 1.009 1.007 .002 
2.0 1.003 1.003 .000 
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Fic. 3. Plot of electron distribution function for 2s electron, D»,=47r°b?, with different types of functions. Li I. 


A, 4-parameter; B, Guillemin-Zener; C, Slater; D, hydrogen-like. 
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Fic. 4. Total electron distribution function for Li I, D=42r*{/*, as a function of r in atomic units. 


electrons, found by using the K functions pre- 
viously described. To calculate this latter quan- 
tity at a distance R from the nucleus, use is made 
of the equation derived from electrostatics for 
the field due to a distribution of electricity of 
density p=a? 


e Z 8x 7* a 
V= ——+— [ redr+ sr ra*dr. (25) 
R R 


0 R 





If ais given by Eq. (21), (25) becomes 
V=—2Z/R+(2/R){1—e-®(1+ER)}. (26) 


As is seen from Table III, there is quite good 
agreement between the field calculated in this 
manner and that of Hargreaves, the deviation 
being about 2 percent at most. This comparison 
refers only to the 1s functions used, and not to 
the 2s function, however. 


EXPRESSIONS FOR N, L, M 


Eq. (15) gives the energy in terms of N, L, M, defined by (13). Algebraic expressions for N, L, M 
for the configuration (1s)? 2s and for the 2s functions of (7) are given below: 


a=a/t; = 5/§ 1’ =n/é. (27) 
2K=f'+9'; b=1+n'; fHl1tsy;. g=1+K. (28) 


3a” a” 


x 


N= |, -576 —+384 ——3 —4+—-—— 
, bs 


n 


72 72 72 72 


. a a a a a 
L=;' 3(2Z — 3) —+(1.5Z-—3) —+3 —+6 —+ (516 —384-Z) ——576-Z — 
| n° n’’ Ke p5 b? 





a a 1 64 
|! (29) 


bafs K4 - fs 
P a” 3840” 


bs d4(2+0)3 
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11520” a’ a’ a’ a’ a’ a’ 
————-+ (4-22) —-—4—-—6 —+(128-2Z—256) —+192Z —+ (384Z — 384) — 
b4(2+0)! K3 2 gt fb bif? bifs 
128a’ 384’ 1536a’ 1536a’ 128a’ 192a’ a’ 
+ + + + + + —3(2Z—§) — 
f(2+b)> f2(2+b)* HF 2+)? oY*2+S/)®> F(f+5)> G+) K* 
+(2Z—§) +(Z- * ams +— ° +-(20- 64Z) — * 4128- 642) - a F ( 
Ds — _- ) 
rae fr f fs fe f(24f)? F2+F) 
of" a” Pa a” n'a 2 a’! tna! eo! a’ ¢! 
M=r' rae PT i i pla —12 en ar a 
n b? bs bs K* K* bf3 bif3 K3 
i a’ n’ 2 (’ 1 | 
+384 ——+-384 ——+-—--+-—- — 125 —— 64 — . (31) 
pp gel 
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